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1 Some properties of materials 

1. The magnitude of things 

Rearranging the things cited, we find that. 

(a) Using E = me 2 for the rest mass of hydrogen, this gives 9 x 10 16 J 

(b) Daily energy consumption of the UK is 3 x 10 15 J 

(c) 1 atom of uranium, when fissioned, produces about 200MeV energy. 1kg of uranium-235, 
when fissioned, produces 200 x 1.6 x 10 -13 x 6.022 x 10 26 /235 = 82 x 10 12 J 

(d) Deuterium fusion releases about 3.65MeV, depending on whether it makes tritium or helium. 
Per kg this is 3.65 x 1.6 x 10“ 13 x 6.022 x 10 26 /4 = 88 x 10 12 J 

(e) The binding energy of H 2 is 4.52eV, and of H is 13.6eV, so to ionise everything requires 
(4.52 + 2 x 13.6) x 6.022 x 10 23 x 1.602 x lO -19 ) = 3 x 10 9 J 

(f) Zero point energy in hydrogen = Tuv /2 and 4161 cm—l=0.51eV, so 50 x 10 6 J. 

(g) Gravitational energy for tidal power mgh = 49 x 10 6 J 

(h) The energy density of coal is 24MJ/kg, so 1kg produces 24 x 10 6 J 

(i) The gravitational potential energy of the sandwich is GMm/r, 12 x 10 6 J. 

(j) A KFC doublicious chicken sandwich contains about 2 x 10 6 J. Source: KFC website. 

(k) The density of air is 1.22%/m 3 , so the KE is 61 x 10 3 J 

(l) Chicken sandwich is a hydrogen-bonded solid, so estimate the specific heat capacity to be 
similar to ice, 2.1kJ/kg/K. If the sandwich weighs 0.2kg then the thermal energy to heat it 
40K is 17 x 10 3 J 

(m) The thermal energy of hydrogen molecules is \RT = 3.7 x 10 3 J 

(n) The kinetic energy of a railway sandwich is \mv 2 = 250J 

So there are some 14 orders of magnitude between the kinetic and rest mass energy of a sandwich 
on a train. No physical quantity has ever been measured with that precision. There are almost 
enough calories in a chicken sandwich to put it into orbit. 

Although the questions mix large and small things, the answers cover a much larger range of 
energies. From this we obtain some idea of the order of magnitude of things. Rest mass is by far 
the largest store of energy, followed by nuclear energy, atomic physics, chemical bonding, gravity 
and finally everyday kinetic energy. Each is separated by about 3 orders of magnitude. 

2. Heating and metabolism 

We will assume that all the energy used in metabolising is converted to heat in the air. 

Amount of heat produced by the students: 

E=nxPxt= 6.48 x 10 7 J 

Assume that this all goes to heating the air 

E = Me AT 

mass of air, volume times density, is M = 1200 kq. heat capacity is mass times specific heat capacity: 
C = 1200 x 1000 = 1.2 x 10 6 J/K 


AT = E/mc = 54 K 
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Final temperature = 20 + 54 = 74°C 

That’s hot! People generate a significant amount of heat, roughly equal to our energy intake, 
around 3000 “nutritional calories” or 12 million Joules per day. 

You might also ask how stuffy it gets: typically people use about 25 litres of oxygen per hour so 
even 100 students make a small impact on the total oxygen content. 

Obviously the sun is much bigger than a person, a more interesting comparison eliminates this 
effect. The mass of the sun is 2 x 10 30 kg, so it produces energy 1.93 x 10 ~ A Wkg~ 1 . By contrast, 
a person is about 60kg, so metabolises lAWkg -1 . Per unit mass, you produce 10,000 times more 
heat than the sun! 

3. Thermal properties in food science Strawberries are 88% water, and water has a high spe¬ 
cific heat capacity so we assume this component will dominate. The specific heat capacity for 
strawberries is then 


Cstrawberries = 0.88 x 4.2 x 10 3 = 369 6J/kg/K 

(a) Energy lost A E = mcAT = 38808 J 

(b) “Respiration” assumes strawberries are emitting heat, so they will need more energy to cool, 
and this will take longer. We don’t have quite enough information, so we need to make a 
sensible assumption about how long the process will take. A fridge uses a few tens of Watts 
(looked up on the internet), so 38808J will take about half an hour assuming reasonable 
efficiency of the fridge(30 x 60 seconds). Doing a sanity check, this sounds reasonable. 

So the heat produced by respiration is 

q = 0.5 x 0.21 x 30 x 60 = 189J 

This is tiny compared with (a), so the effect can be safely neglected. Errors in the half-hour 
assumption do not change this conclusion. 

(c) Strawberries are mainly available in summer, and since this question comes from Edinburgh 
we’ll take 25°C to be room temperature. From (a) we need to produce 38808J to heat back 
up, and since they are thermally insulated the only source is from respiration which takes: 

time = 38808/0.5/0.21 = 3.7 x 10 5 s 

which is over four days, consistent with our conclusion in (b) that respiration is generally a 
minor effect compared with heat flow. 

(d) Assume the respiration (105mW) is related to the breakdown of the strawberries, and loss 
of nutritional value (700kJ). This amount of energy is given of in heat in (700 x 10 3 /0.105) 
seconds, which is about 83 days. Assuming that 10% turned-to-mush is unacceptable for 
human consumption, this gives a lifetime of just over a week. 

4. Phase Change: latent heat 

To do this question we will make a strong assumption: that we can reach the answer by considering 
only the conservation of energy between initial and final states. This is justified because the 
insulated container allows no heat to enter of leave (and no work to be done). 

a) We must consider whether cooling the water 0 °C would release sufficient energy to provide the 
latent heat of melting the ice. 

The latent heat required to melt the ice is given by 

A E l = 334 x 0.03 = 10.02fcJ 
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The available heat in the warm water as it cools to 0 °C is 

A E = 4.2 x 0.2 x (20 - 0) = 16.8/cJ 


So there is more heat available than needed to melt the ice, and all the ice will melt. There is 
6.78kJ of excess heat remaining, so the final temperature is 


T = 0 + 


16.8 - 10.02 
4.2 x (0.2 + 0.03) 


= 7 °C 


Notice that if A E — A E L < 0 then this equation gives a water temperature below zero. This is 
unphysical: the criterion for the equation to be valid is that all the ice melts. We could blindly use 
the equation, then if the temperature turns out to be positive it proves the assumption was OK. 

Also notice that the mathematical steps do not respect any physical process. They imply first 
cooling the original water to 0°C, and putting the energy “somewhere”, then taking some of the 
energy to melt the ice, and finally using the remaining energy to heat all the water. 

The system has moved towards equilibrium. There is nothing in the calculation which tells us this: 
we used our physical intuition to say that heat should flow from the hotter body to the cooler 
body and the final state should be at the same temperature throughout. Later, we can replace this 
intuition with the Second Law of Thermodynamics. 

5. Gravity and heat 

Assume that all the gravitational energy is converted into heat energy. 


mgh = mcAT 


Thus 


AT = gh/c = 


9.8 x 50 
4200 


0 . 12 # 


Of course, not all the gravitational energy is turned into heat. Sound is produced. Erosion occurs 
by breaking of chemical bonds. Spray and falling water exchange heat with the atmosphere and 
absorb heat from the sun. 

Although most of the water in the Niagara River is diverted from the falls to generate hydro electric 
power, this does not change the conclusion regarding temperature: the volume of water does not 
appear in the calculation. 

In 1847 Joule attempted to demonstrate the mechanical equivalent of heat by measuring the tem¬ 
perature differences in a waterfall (Cascade de Sallanches, near Chamonix). He failed. 

6. The ideal gas law 

Using a 0 subscript to define that standard Pq^Tq are 273K, 101325Pa, for one mole of ideal gas 
we have Vo = 2.2414 x 10 _2 m 3 , n=l. 

Our actual sample has T = 298#, n = 1, P = Po and, 

D _ P 0 Uo _ PV 
nT 


2.447 x 10 ~ 2 m 3 

The ideal gas expands when heated. 


nT 0 


whence 


Vb 

V 


n 

T 


V = —2.2414 x 10“ 2 = 

973 
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The ideal gas features highly in thermodynamics, and surprisingly turns out to be a good model for 
many things, including dilute solutions and cavity radiation. But it fails at low temperature, due to 
interatomic forces and quantum mechanics. For that reason, chemists use “standard temperature 
and pressure” as a reference state, rather than “zero”. 


7. Molecular ideal gas 

For an ideal gas PV=nRT, where n is the number of molecules , not atoms. So 


PV 1.2 x 101325 x 82 x 1(T 6 


RT 


8.3 x 300 


= 0.004 mol 


O 2 is a diatomic molecule, and has a molecular weight of 32, which is the number of grammes per 
mole, so converting to SI units. 

m = 0.004 x 32 x 10 -3 = 1.28 x 10 ~ A kg 

For O 3 , the molecular weight is 48, so the sample mass would have been 1.92 x 10 ~ A kg. Note that 
if the ozone decomposed into O 2 (at constant volume and temperature) then the pressure would 
increase by 50%. 

In classical physics, every atom at equilibrium has the same amount of energy. This is called 
“equipartition”. Some of this goes into molecular vibrations, which leaves less available for mov¬ 
ing the whole molecule to bang against the contaner walls, the microscopic process which causes 
pressure. This is why a gas of atoms arranged as O 3 molecules has lower pressure than teh same 
atoms as O 2 - more of the energy is in bonds which do not contribute to the pressure, (n.b. in a 
solid the pressure does include the energy needed to compress the bonds.) 

8 . Melting, heating and boiling 

The heat is supplied at a constant rate of lkW, so the time taken in seconds will be equal to the 
number of kJ of energy required, so, 

a) Heat ice: mcpAT = 1 x 1.94 x 4 => 7.8 seconds 

b) Melt ice: ra/ me /t = 1 x 334 => 334 seconds 

c) Heat water mcpAT = 1 x 4.2 x 100 => 420 seconds 

d) Boil water ml vap = 1 x 2256 => 2256 seconds 

e) Heat steam mcpAT = 1 x 10 x 2.04 => 20.4 seconds 

Melting involves structural rearrangements, breaking some bonds to change a rigid structure to a 
fluid one. 

Boiling requires much more energy as all bonds between molecules must be broken. 

Heating adds thermal energy to individual molecules, including vibration in ice, motion of molecules 
in water and steam, and vibrations of the atoms in the molecule itself. 

Comparing boiling with melting, we see that about l/ 8 t/l of the bonds are broken in melting, and 
this is enough the allow the water molecules to move about over time. In liquid water, bonds are 
continually breaking and reforming, sometimes between different molecules. 

9. Conduction of Heat 

Strictly, this question is outside the realms of thermodynamics, as the heat conduction equation 
cannot be derived from the Laws of Thermodynamics. 

a) Consider the system at time t where T(t ) < To. In a small interval of time At, the heat transfer 
into the system is 

io= !^skfim At 

New temperature is then 

T(t + At) = T(t) + 
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so that 

KA(Tq — T(t)) . 
me AT =-- y -^At 

T 

assumptions being that AT << T, heat is only lost through area A, the heat source has a homo¬ 
geneous temperature throughout. 

b) In the limit of small At, can convert to derivatives... 


dT _ KA[T 0 - T(t )] 
dt mcL 


l 


dT 


In 


i T-T 0 
T\ - T 0 


KA 

mcL 

KA 


/ 


dt 


mcL 


(h ~ h) 


T 2 -Ti = (T 0 -Ti) 


1 — exp ( — 


K Afo-h) 
mcL 


Sanity check: as oo, T<i To. eventually, the object reaches the temperature of the heat 

source reservoir. 

c) Aluminium 39.53°C, 60°C; 

Porcelain 35.04°C 54.4°C; 

Rubber 35.004°C 38.48°C; 


Wool 35.001°C; 36.2°C. 

Aluminium insulation is a bad idea. 
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2 Temperature scales, work, equations of state 

1. Temperature scales: influence of thermal properties 

The scales will not necessarily give the same answers. 

In order for alcohol and mercury thermometers to agree at all points the thermal expansion must 
change in the same way with temperature, because the “thermometer” reading is related to the 
thermal expansivity of the material. If the thermal expansions were constant with temperature 
then if they match at of 0 °C and 100 °C they will agree with each other and with the Celsius scale 
at all temperatures. The actual value of the thermal expansion does not matter. 

If thermal expansion varies linearly with temperature, they will agree with one another, but not 
with the Celsius scale. 

If the thermal expansions varies nonlinearly with temperature, the thermometers will not agree, 
unless there is some freakish cancellation of errors. 

Note that the thermal expansion applies to both the mercury/alcohol and the materials making 
up the thermometer. And further if the pressure of the mercury/alcohol changes then the bulk 
moduli must also be considered. 

2. Temperature scales: based on electrical resistance 

Roger’s “Resistance temperature” (Tr) is defined relative to the ice temperature 273.15K 

Tr = kR = kR 0 ( 1 + aT + /3T 2 ) 

where T = T^eaZ ~ 273.15Ah This guarantees that the scales match at the ice temperature. We 
need to find the constant k so that the scales match at the triple point (T = 0.01): 

= Tj d l al = kRo(l + a x 10 “ 2 + (3 x 1CT 4 ) = 273.16# 


273.16 

i? 0 (l + a x 10- 2 + 0x 10- 4 ) 


Now use this value of k 


r C = kR{™°C) = TY- 1 


(1 + a x 70 + /3 x 70 2 ) 


+ a x 10 “ 2 +(3 x 10 -4 
A couple of lines of code will give you the discrepancy at all temperatures. 


341.79# 
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The maths is easy, but the physics is quite subtle. Practical thermometers must use some physical 
object to measure temperature. The notion that nature has some absolute scale for temperature 
is the Zeroeth Law of Thermodynamics. 

Linguistically, the plural of Kelvin is Kelvins, while for Celsius it is degrees Celsius. 

3. Work Done in various processes 

(a) Ice has lower density than water at CPC, so the volume decreases on melting. Assume pressure 
to be constant, Work = -PdV 

rwater 

W = - / PdV = -P a t m x(M/p water -M/p ice ) = -1.01 xl0 5 x[(10/916) - (10/1000)] = 92 

J ice 

This is positive, which means work is done on the system - the pressure favours melting. 

If the volume of a system expands, it has to do work on the atmosphere to push it back. Here 
the opposite is true, the atmosphere presses on the system and does work as it contracts. 
Strictly, latent heat depends on whether we have constant pressure or constant volume melting. 
They differ by the work done, as calculated above. However, this energy is tiny compared to 
the latent heat (3340 kJ) so the two latent heats are essentially the same. 

(b) Make the analogy to obtain an expression for work done in stretching, dW = FdL. Check 
dimensions: Force x length = Energy. Consider the sign: doing work on a gas would compress 
it, hence —PdV doing work on a wire would extend it, hence -\-FdL. 

W = J FdL = F(L 2 — Li) 

The tension in the wire in F=mg=20N, and the length change is 0.0001m, so the work is 
0.002J. 

For the thermal expansion, we write the temperature change in terms of two variables: 


dT = 




dF 


At constant force, dF = 0, and we are given that the linear thermal expansion is: 


0 = 


1 

L 



= 16.6 x 10" 6 


So we can write: 


P-Jn 

AT = ln(1.001)/16.6 x 10 -6 = 60K 

The required temperature change is fairly substantial 0.001/16.6 x 10 -6 =60K. 

n.b. given that the length change is small, you can get the same answer using: ln(^) ~ 

The hidden assumption here, that linear thermal expansion coefficient is at constant force (or 
tension). Does the mass matter? 

The heat capacity of copper is 390J/kg/K, density 8960 kg/m 3 . Supposing the wire cross 
section is 1 mm 2 , then the mass is 9 x 10 -4 and the heat supplied is 21 J, far in excess of the 
work done lifting the mass. So in these conditions the linear thermal expansion coefficient is 
independent of the tension. 

(c) The electrical work is V 2 t/R =(36 x 10)/(4.2 x 10 -3 ), a fairly hefty 86kJ. Enough to melt the 
wire. What we are doing here is, essentially, creating a short-circuit and blowing a fuse. 
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This question illustrates that mechanical work in condensed phases is typically small in comparison 
to the thermal energy at room temperature. So much so that reference quantities like “latent heat”, 
“Young’s modulus”, “resistivity”, even “heat capacity” for solids are usually quoted without stating 
whether they are constant volume or constant pressure. 

4. Calculating Properties from the equation of state, and vice versa 

(a) First, we show that measurable material properties can be calculated if we know the equation 
of state. 

For properties of the ideal gas, we proceed by differentiating the equation of state PV=nRT. 
Definitions: isothermal compressibility k = ^ (§rp) T ; isobaric expansivity /3 = y (f^) P 

Notice that these are log derivatives k = — ( d gp / ) T i P = Ultimately this is because 

we are interested in fractional change in volume, not absolute change. 

We rearrange the equation of state to make V = nRT/P and evaluate the differentials, then 
use the equation of state again to simplify the final expression. 

1 /<9U\ _ 1 nR _ 1 

V [dfjp ~ V^P ~ T 


_ -1 (dV 
V (dP 


T 



1 

P 


By definition, the bulk modulus is the reciprocal of the compressibility, so B = 1/k = P. 

The question illustrates the central role of the equation of state for determining material 
properties. The myriad of measurable physical properties are not independent, they can all 
be calculated from the equation of state. 
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(b) Now we don’t know the equation of state: we must work it out from materials properties. So 
we reverse the process in the previous section, integrating instead of differentiating. So for 
this substance, 


<9U\ 

dp) T ~ a 
J dV T = J -adP T = V = -aP + f(T) 

Where f(T) is some as-yet unknown function of T only. 


( d T) 

\dT J j 


2 bT 


J dV T = J 2bdP T = V = bT 2 + g(P) 

Where g(P) is some as-yet unknown function of P only. 

Now we match the two results to obtain 

bT 2 + g(P) = —aP + f(T) 

bT 2 — /(T) = —aP — g(P) = constant 

This must be constant, since T and P are independent. So 

V = bT 2 — aP + constant 

Notice that knowing the properties doesn’t fully specify the equation of state. In this case, 
the undetermined constant of integration is the volume at zero temperature and pressure. 

(c) Finally, we show that some materials properties can be calculated directly from others. 
Remember Cp — Cy — R and PV = RT for one mole of an ideal gas. 

The volume of a given system containing n moles of a monatomic ideal gas is given by V = 
V(P,T ), where P is the pressure and T the temperature, and PV = nRT. 

_ PV _ V/3 _ VT(d 2 
T K K 

In fact, this is a general result for any material. The difference in heat capacities is the extra 
work needed to expand the material in the constant pressure process. This work depends on 
how much expansion we get, /?, and how resistant the materials is to expansion. 

This illustrates that thermodynamic properties are not independent. It is sometimes possible 
to work out a hard-to-measure quantity by measuring other material properties. 

Notice also that the compressibility and expansivity are independent of the amount of sub¬ 
stance. 

5. Temperature scales: based on water 

Using p = A + BT + CT 2 + DT 3 . 

We are trying to fit the isobar for water to a cubic equation. We could directly set up four equations 
in four unknowns, but there is a simpler way to proceed. The equation is quartic, so we can define 
the zero of temperature for mathematical convenience. Once the parameters are determined, their 
values for other choices of the temperature are just linearly shifted. 
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If we define a temperature scale where T=0 equates to the maximum density, then it follows 
immediately that A= 1000 kg/m 3 and B = 0 

Now the low temperature data point is at T=-44, the high temperature one is at T=96, so we have 
two simultaneous equations: 


960 = 1000 + 44 2 C - 44 3 L> 

960 = 1000 + 96 2 C + 96 3 D 

From which we obtain: C=-0.0125 D=0.000044. this is true for the temperature scale with zero at 
+4°C. If we wanted to use T in Kelvins then the equation of state would still be quartic: 

p = A + B(T — 277) + C(T — 277 ) 2 + D(T — 277) 3 , but the power series in TKelvin would have 
coefficients Ax = A — 277B + 277 2 C — 277 3 D etc. 

for the present calculation, retain the scale with zero at the maximum density. 

In the thermometer, the water column rises by 1 %, but the diameter is unchanged. So the volume 
expansion (or density decrease) required is 1%, to 990 kg/m 3 

-10 = —0.0125T 2 + 0.00044T 3 => T = -27, +30 

Converting back to Celsius, this gives two solutions, -23°C or +34°C. A sketch of the isobar will 
confirm this is reasonable. 

Even assuming the water can be prevented from freezing, density of water is a terrible choice 
for a thermometric property because a single value for the density does not uniquely define the 
temperature. 

6. Joule’s experiment 

(a) Assume that the Gravitational Potential Energy ( mgh = 20 . 10.2 = 400J) is converted to heat: 
mcAT = mgh. 

AT = gh/c = 10.2/4200 = 0.047 K 

(b) Some energy is still in kinetic. KE = \mv 2 = 2Qx 2 Q ‘ l2 = 0.1J. But only one part in 4000, or 
0.025% 

This calculation follows the demonstration by Joule that potential energy could be converted to 
work, and then to heat. In the actual experiment great care has to be taken to avoid or compensate 
for heat loss, and there was some scepticism about Joule’s claim to be able to measure temperature 
to the required accuracy. 
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3 Work and heat, the First Law 

1. Work done in the expansion of an ideal gas 

In an expansion, dV > 0 , so positive work is done by the system, negative work is done on the 
system. It is best to think about the sign of the work, rather than try to memorise equations with 
or without minus signs. 

The internal energy of an ideal gas depends only on the temperature, so it can’t change in any 
isothermal process. As a consequence, in parts (a) and (b), if work is done, and equivalent amount 
of heat must be added. 


(a) Start by writing the general expression for volume: 


dV (P, T) 



fdV 

\dP 


dP 

T 


notice that the first term is zero for isothermal process (dT = 0). 
equation 


dV = 



dP = - 

T 


nRT 

-p2~ 


dP 


then using the ideal gas 


So that 


dW = -PdV = 


nRT 

P 


dP 


(b) Gas expands isothermally (T = T 0 ) from initial pressure Po- PoVo = tiRTq = nRT = PV. 
since V = 2Vo, then P = P 0 / 2 . Now using the result from (a) 


/ r m 1 

PdV = / -p-dP = nRT [\nP]p° o /2 = nRT In - = -nRT In 2 

Notice that there are three minus signs here - one from the definition of work, one from 
differention, and one from the logarithm. 

(c) Following the same logic, this time for isobaric process (dP = 0) 


dv ^ = {w) p dT ='£ dT 

W = - J PdV = -nR J dT= -nRAT 

You need to be careful about the definition of work. There is nothing in nature to define 
its sign, indeed ehemists and engineers tend to use opposite definitions because ehemists are 
interested in what happens to the material, whereas engineers are interested in the work the 
engine does. 

2 . Work, heat, PV diagrams 

Every material has an equation of state, which can be written as T(P,V). Thus a point on a 
PV diagram is sufficient to give all the information about a particular sample of known mass at 
equilibrium. Any process is then fully represented by a line on a PV diagram. We could instead 
use PT or VT diagrams, but PV proves most convenient because of the direct link to work: f PdV. 

(a) The First Law tells us that the change in internal energy comprises the heat put into the system 
(+80J) minus the work done by the system: 80 — 30 = 50J = AU = Ub — Ua 

(b) U is a state variable, so its change is independent of the path taken: A Uacb = A Uadb = 
A Uab — 50J. Since the work done on the ADB path is only Wadb — 10J, the balancing heat must 
be 60J. 
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(c) Again, for a state variable on the reverse stroke, A Uba = —A Uab = — 50 J. Now the work is 
done on the system, the sign is also changed, so the heat going into the system must be —50 — 20 = 
—70J. Of course, negative heat going into the system means positive heat given out. 

(d) A Uad = 40J, so AUbd = A Uab ~ A Uad 10*7 


For the work: Wadb = Wad + Wdb — —10J. but since BD is a constant volume process, no 
work can be done, so Wdb = 0 and Wad = —10 J From which the heat inputs can be seen to be 
divided Qad = 50 J, Qdb — 10 J. 

The question illustrates that there are different ways to get from one state to another, involving 
different amounts of work and heat input. However, for a given start and end point, the sum of 
work and heat is independent of the process. This is a concept on which many thermodynamics 
calculations are based - to determine the properties of the final state of a sample, we can consider 
any path to get there, including irreversible one. Samples at thermodynamic equilibrium have 
no memory. The inverse is also true, anything with a memory connot be at thermodynamic 
equilibrium. 

3. Free expansion of van der Waals gas 

In a free expansion, no work is done on the surroundings, and no heat flows into the system. 
Consequently, according to the First Law, the internal energy U must remain constant. 

Given that 

fdT\ _ _ a /n\ 2 

\dv) u ~ “cv vyJ 

which represents the free expansion coefficient for T changing with V. We can get the temperature 
by integrating from the initial state (1) to the final state (2)... 



T 2 


Ti = AT = - 


2 

an 

Cy 



v 2 

Vi 


an 2 / 1 1 \ 

c^\w~v 1 J 


We know Vi, but it’s not easy to extract V 2 . But we do know that a gas (state 2) is a lot less dense 
than a liquid (state 1), so we will assume that it is big enough that - W 

Now, based on this assumption, The initial pressure comes from the van der Waals equation: 

n\RT\ 


Pi = 


Vi — nib 


- ^ = 2.49 x 10 7 Nm~ 2 


compared with one atmosphere (P — 2 = 1.01xl0 5 Nm 2 ) we see that Pi >> P 2 . For a gas, this 
implies that our approximation to neglect ^ will cause less than 1% error. So finally 


AT 


an 2 

CvVi 


-2.29 K 


Recall that the constant a represents binding between atoms in the VdW gas. When you expand 
a VdW gas, this binding gets weaker. That energy has to come from somewhere, and in a free 
expansion the only possible source is the thermal energy. 

4. Free expansion experiments and internal energy of ideal gas 

The observation can be written mathematically as: 



u 


= 0 
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we can use the reciprocity relation 


/ dT\ fdU\ fdV\ 

{WJ u {df) v {m) T 


-l 


to give. 



so that the experiments imply that one of those partial derivatives must be zero for this gas. 

Now (§^) v = C v , and we know that heat capacities must be finite, so the measurement implies 
that... 



i.e. if we expand the gas at constant T, the internal energy doesn’t increase: U is independent of 
volume U = U(T) only. Remembering the model of the ideal gas as hard, elastic particles, it is 
reasonable that their only energy is kinetic energy, i.e. temperature. 


dU = 


(dU 

V<9R 


T iv+ {%) v iT - cviT 


5. Adiabatic processes 

For an adiabatic process, dQ = 0 so the first law is dU = —PdV , or using the ideal gas specific heat 
C v dT = —(PdV). The differential form for n moles of the ideal gas equation gives PdV + VdP = 
nRdT. We are looking for the relationship between P and V, so we combine the equations to 
eliminate dT\ 


dT = 


PdV + VdP 
R 


P 

TlCy 


dV 


Rearranging: 


dP _ nRpC v dV _ dR 

~p ~ cl V " _7 V“ 

Where we used the fact that nR = Cp — Cy and defining 7 = Cp/C v . 

(n.b. you can do this by considering one mole of ideal gas, in which case n = 1 and C v = cy) 
Thus ln(PV 7 ) is constant, and so is PV 1 . 

7 is a ratio of specific heats, so it has no units. The units of c depend on 7 . e.g. for monatomic 
ideal gas 7 = 5/3 and c has units Jm 2 , for a diatomic gas it is Jm 4 . Notice that although 7 is 
fractional, the units always have integer dimensions. 

It is very important to define a sign convention with respect to work done. 

W = -PdV as the work done on the gas. (think: positive work means negative dV, energy needed 
to push a piston to compress the gas) 

W = +PdV as the work done by the gas. (think: positive work means positive dV, energy from 
the gas pushes a piston out) 

Here we are interested in work done by the gas, so 

r v2 

W= PdV 
Jv 1 
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For adiabatic processes, we know that PV 7 is constant. Calling the constant “c” and eliminating 


P. 


Wmc 


rV2 

/ 

Jv i 


V~ 7 dV = c 


Y 1-7 1 V2 


1-7 


Vi 


1-7 


(V l-7_ V l-7) 


now noting that c = P\V\ — P 2 V 2 we obtain 


PiV 1 -P 2 V 2 


Sanity checking, if P\V\ = P 2 V 2 (i.e. T\ = T 2 ) then no work is done, which seems sensible since 
internal energy U(T) is unchanged and adiabatic means no heat is exchanged either. 

If 7 = 1 infinite work would be done. Of course, for an ideal gas 7 = c p /cy is never 1 , but it gets 
closer for a polyatomic ideal gas. Broadly what is happening in a polyatomic-molecular ideal gas 
is that much of the energy is tied up in molecular vibrations, and a smaller fraction in translations. 
Only the translations contribute to pressure, as the molecules bounce off the container walls. So 
for a given change in temperature, extra energy is released from molecular modes to do work. 

6. Isobaric processes 

From the First Law, A U = AW + A Q 


AQ = C p AT = C p (T 2 -T 1 ) 

AW = -PAV = -P 0 AV = -P 0 (V 2 - Vi) 

so that 

AU = U 2 -U 1 = C p (T 2 - Ti) - P 0 (V 2 - V{) 
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This is completely general, now we consider the special case of an ideal gas. We can use the 
equation of state PV = nRT to eliminate temperature 


U 2 -Ui = %{P 0 v 2 

nR 


PqV\) - P 0 (V 2 


«>-<S 


- 1)P 0 (V 2 - Vi) 


We could plug and chug the numbers at this point, but we can go a bit further in the algebra with 
Cp — C v = nR and 7 = Cp/Cy. So that 


Cp 

nR 


r£-i = (7-i) 


-1 


so that 


U 2 -U 1 = Po ^ 2 TA = l x 1.01 x 10 5 x (10 - 5) = 7.6 x 10 5 J 
7 — 1 2 

Notice once again the pivotal role of 7 — 1, which microscopically depends on the number of degrees 
of freedom per molecule 7 = (m + 2 )/m. 

Also notice that the heat in an isobaric process goes to internal energy and work A Q = A (U + PV). 
It can be convenient to define a state variable H = U + PV whose changes are equal to the heat 
supplied. This state variable is Enthalpy, and is used a lot in thermodynamic descriptions of 
isobaric processes. This is a recurrent theme in thermodynamics, although we can use any state 
variables to describe the system, choosing one which matches the boundary conditions makes things 
much easier. 

7. Cooling in an adiabatic expansion 

To derive the Joule-Kelvin expression, start with the triple product 

fdT\ __(dT\ ( dH \ 

\dP) H ~ \dH) P \dPj T 

We can now identify the specific heat at constant pressure Cp = (§^) p , and the given expression 
for (§^) T - Substituting these into the triple product gives: 



H 


1 

Cp 



1 

Cp 


(TV0 - V) 


where (3 is the thermal expansivity. You can check your working by noting that (§p) H is zero for 
an ideal gas. 

To use this for the Joule-Kelvin cooling, we can always write T in terms of any two state variables, 
e.g. T(P, iJ), from which 


dT = 



dP + 

H 



dH 
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By definition, for an isenthalpic process, we have dH = 0, so 


dT = 


fdT 

\dP 


H 


dP 


and provided the Joule-Kelvin coefficient (§p) H is independent of T we can separate variables and 
integrate. 


A T = 



dP 


We can anticipate that the adiabatic expansion produces more cooling than a throttling process. 
If we look at the PV diagram for an ideal gas, we see that adiabats P oc V~ 7 drop more steeply 
than isotherms P oc E _1 , since 7 > 1. This means that in an adiabatic expansion we drop from 
a higher temperature to a lower one. So if V\ < V 2 then T\ > X2, the gas cools as it expands 
adiabatically (no heat is exchanged, but work is done by the expanding gas, and the energy to do 
this work must come from lowering the internal energy). 

For an ideal gas, internal energy depends on T only, so a Joule process (constant U ) follows an 
isotherm. The definition of Enthalpy is H = V + PE, and for an ideal gas the product PV is 
a function of temperature only, i.e. PV = nRT. Thus for an ideal gas the isenthalp is also an 
isotherm. 

For a material “similar” to the ideal gas, we can assume that the isenthalp remains close to an 
isotherm. 



To prove the inequality for any substance, consider the triple product for the relevant quantity in 
an adiabatic process: 


Apply a Maxwell relation 




p 





p 


Applying the definition of C pj the second two terms become: 


(dT_ 

\<9P 


s 


fdV\ T_ 

\df) P cp 
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To see the difference in cooling between the isenthalpic and adiabatic expansion process, we look 
at: 


so that 


A Ta - A T jk 




dP 


A Ta - A Tjk 



" T fdV\ 
Cp \df) P 



t rov\ 
Cp \df) P 


v 

Cp 


dP 


and since P 2 < Pi in an expansion, while V and Cp are always positive, this integral must be 
negative. Both A Ta and A Tjk are negative on cooling, so the isenthalpic (Joule-Kelvin) throttling 
process produces a smaller cooling effect than the adiabatic process. 

A final note: William Thomson became Lord Kelvin in 1892 (strictly, Baron Kelvin of Largs). 
Some people still haven’t noticed and refer to isenthalpic expansion as the Joule-Thomson effect. 
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4 Cycles and the Second Law 

1. Statements of the Second Law of Thermodynamics 

Clausius : It is impossible to construct a device that, operating in a cycle, produces no effect 
other than the transfer of heat from a colder to a hotter body. 

Kelvin-Planck : It is impossible to construct a device that, operating in a cycle, produces no 
effect other than the extraction of heat from a single body at a uniform temperature and performs 
an equivalent amount of work. 

Suppose that the Clausius statement of the Second Law of Thermodynamics is false: heat can flow 
spontaneously from a cold body to a hot body. Label the heat transferred by the Clausius-violating 
device Q i. Now connect a separate standard heat engine to the two reservoirs, and adjust it so 
that it takes heat Qi from the hot body (dumping heat Q 2 to the cold body, and producing work 
W = Qi — < 32 )- Not net heat is taken from the hot body, so taken together, the Clausius-violating 
device + standard heat engine transforms heat Q\ — Q 2 from the cold reservoir directly into work. 
This violates the Kelvin-Planck statement of the Second Law. 

This may also be proved by adjusting the standard engine to dump heat Qi to the cold reservoir: 
now the heat come from the hot reservoir. 

Notice that all the devices considered here obey the First Law of thermodynamics. 


2. Efficiency of engines part 1 


Efficiency of a heat engine, by definition, is 77 = 

W 


T] = 


Q1 


(work out)/(heat in). Here 



For the special case of a Carnot engine 


Vc = 1 - 


T2 

Ti 


To find which increases the efficiency more, we need to compare the change of increasing the hot 
temperature dr\ c jdT\ with reducing the cold temperature —dr] c /dT 2 . Notice that both of these 
quantities are positive. 


drjc 

dT 1/T2 


Since, by definition, T-\ > T 2 , it follows that (—fjjjj 


dr)c\ 

, _1 


dT 2 ) 

Ti T, 


t 2 

( 0r h ' 

\ x t 2 

T\ 

A dT 2 , 

> Tl ?! 

that ( 

drjc \ 
dT2 ) Tl 

> ($?) 


t 2 


so a decrease in the cold 


reservoir is more effective in increasing efficiency than an equivalent rise in the hot reservoir. 


3. Efficiency of engines part 2 

For the given system to maintain the cold reservoir at X \ no net heat should enter or leave it, hence 
the two Q 2 flows with opposite sign. 
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From the first law, the work coming from the engine W\ = Qi — Q21 while the work going into 
running the fridge is W a = Q a ~ Qz- For a compound device where the engine drives the fridge, 
the net work out is 


W = W 1 -W a = (Ql- Q 2 ) - (Qa - Qi) =Ql-Qa 

So the efficiency of the compound engine is work done divided by heat input. 

= w = Ql - Qa 
V Ql Ql 

at best, with reversible Carnot engine and fridge, this would be, 



the same as if we ran the engine between the hot reservoir and the ambient one. 

So maintaining the cold reservoir doesn’t gain us anything. You can’t beat the second law. 

4. Efficiency of engines part 3: The ’Otto cycle’ 

Process ab is an adiabatic compression of the air as the piston moves 
up, work being done on the system. 

Process be is a const ant-volume heat transfer to the air from an ex¬ 
ternal source while the piston is at the top. This process is represents 
the ignition and explosion of the fuel-air mixture. 

Process cd is an adiabatic expansion, the power stroke as the piston 
moves down doing work. 

Process da completes the cycle by a const ant-volume process which 
represents the exhaust. 

By inspection T d > T a and T\ b < T c (convince yourself this is true) 

For analysis, let’s imagine everything is reversible, the working substance is one mole of an ideal 
gas, with internal energy cRT ( c = 5/2 for a diatomic gas, like air). Assume also it is not expelled 
and the heat for the “explosion” is actually supplied externally. 

ab, work is done is on the system, no heat is exchanged, so work is equal to the change in internal 
energy 

W ab = cR(T b — T a ) 

be: heat is added from some external source and the working substance pressure increases at 
constant volume. 

Q bc = eR(T c — T b ) 

cd: Work is done by the system as it expands, again equal to the internal energy change. 

W cd = eR(T c - T d ) 

da: heat is dumped from the system at low temperature 



Qda cR(T d - T a ) 

Net work done by the system is eR(T c — T d ) — cR(T b — T a ) = cR(T a + T c — T b — T d ) 

Note that this is per mole of gas - obviously the actual work done is propotional to the volume of 
the cylinder. 

Remembering that the work done is given by f PdV we can write the work done in the cycle as. 

pb pc nd pa 

/ PdV + / PdV + / PdV + / PdV 

J a J b J c J d 
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Deleting the be and da segments for which dV is zero, and integrating ab in the opposite direct 
gives: 


rd na 

/ PdV- / PdV 

Jc Jd 

Remembering than an integral is just the area under a curve, we can recognise this quantity as the 
area enclosed by the curve representing the cycle in the PV diagram 

Efficiency is defined as the ratio of heat input to work output. For the Otto cycle the work done 
is in the power stroke expansion (cd), less the work needed to recompress the gas (ab) W cd + W a b- 
The heat input occurs in the (be) section. 

Thus being careful about the signs of the work done on and by the system, the efficiency is 

work.done _ cR(T c - T b ) - cR(T d - T a ) T d - T a 

heat.in cR(T c — Tb) T c — Tb 

This is less than for a Carnot engine operating between T c and T a . 

5. Efficiency of engines part 4 

The venture capitalists should consult a competent physicist. Who would check the laws of ther¬ 
modynamics. 

The first law requires conserving energy round a cycle: dQ+dW =0, so 

5000 J = 3500 + 1500 J 
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So far so good... 


The second law implies that the efficiency of the engine cannot be higher than the Carnot efficiency. 
The claimed efficiency is... 


Qout _ i 3500 

~Q~ = ~ 5000 


= 0.3 


While the Carnot efficiency for an engine running between those reservoirs is: 


rj c = 1 - 


^ cold 
Thot 


= 1 - 


300 

400 


= 0.25 


So it is clear that this device is claimed to violate the laws of physics. The physicist will advise 
against investment. 

And the venture capitalists should try to find a “greater fool” and sell on the business 

6. Efficiency of Fridges 

The “efficiency” of a fridge is defined by (heat removed from cold box)/(work in), 
so for a Carnot fridge we have: 


Q 2 = T 2 
Qi -Q 2 Ti - t 2 


If we take room temperature to be 293K, and the cool box 5K above freezing at at 278K, then a 
Carnot fridge should have 


278 

293 - 278 


18.5 


This is usually referred to as the “coefficient of performance”, because an efficiency above 1 doesn’t 
sound right. Note that cooling the ice box is a lot less efficient. Real fridges tend to have v ~ 5. 


7. Efficiency of heat pumps 

The electric heater simply supplies 20kW (i.e. 20kJ/s). 

A Carnot heat pump operates with coefficient of performance Qi n /(Qm — Qout) = Tin/ ( Ti n — T out ). 
In this case it would be 300/20=15, so in principle the heat pump could supply 15 x 20 = 300kW. 
In practice a coefficient of performance of around 4 is typical. 

The heat pump really can outperform an electric heater. Although the electrical work can be 
converted 100% into heat, it is also possible to exploit the temperature difference to generate more 
energy. 

8 . Multipurpose device. 

This device is acting as both a refrigerator and a heat pump, essentially putting the “waste heat” 
from the fridge to good use. Start by estimating the Carnot efficiency of such a device: a sensible 
assumption is needed for the temperatures, so lets take the temperatures as 275K (fridge) and 
295K (living room). These are slightly different: 


r)™ ax = 275/(295 - 275) = 13.5 
r= 295/(295 - 275) = 14.5 
Now look at the claimed efficiencies. For the refrigerator: 


r]R = Qc/W = 400/100 = 4 
Vhp = Qh/W = 1000/100 = 10 
So we see that the claimed efficiencies are perfectly plausible. 
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There is an oddity here, since W + Qc Qh • Does this mean that the device breaks the first law 

and the claim is invalid? 

It doesn’t. In any device there will be losses, and in this case the fridge isn’t anywhere near 
the Carnot efficiency: probably heat is being extracted from the kitchen environment by cold 
pipes outside the fridge compartment. If the company had claimed tjr = 9, extracting 900 W 
of heat from the fridge compartment, then all the energy would be accounted for. Since this is 
thermodynamically plausible, the more modest claim of t]r = 4 certainly must be. 

Don’t worry if you got this wrong, almost everyone does! Coefficient of performance is a tricky 
concept: A fridge where the hot and cold baths had equal temperature would have infinite coefficient 
of performance. All that means is that no work is required to move heat across. If the “hot” 
bath has lower temperature than the “cold” one, the calculation gives a negative coefficient of 
performance. What that means is that heat will flow spontaneously in the direction required. If 
you open a window when its hot inside, then you are using an “air conditioner” with better than 
infinite coefficient of performance. 

9. Yet another cycle 

Let there be n moles of gas in the system, such that the ideal gas law is PV=nRT. 


The system has three states: 

(a) V a , P a , T a , internal energy U = nRT a = P a V a J 

(b) V b ,P a ,T b , internal energy U = nRT b = P a V b 

(c) V&,P6,T c , internal energy U = nRT c = P b V b 


We will consider a sign convention using work done on and heat supplied to the working fluid of 
the system system. 

(i) Work (positive) is done on the working fluid while compressing the system W{ = P a (V a — V b ). 
The heat absorbed this constant pressure process is, by definition, Qi = C p (T b — T a ). This is 
negative - so the waste heat given up is Qi = C p (T a — T b ) 

(ii) No work is done in this isovolumetric process, thus Qa = dU bc = C v (T c — T b ). This is the heat 
input to the cycle 

(iii) No heat is exchanged in this adiabatic process, thus Wm = dU ca = C v (T a — T c ) 

The efficiency (work done)/(heat input) can be conveniently written as 1 - (waste heat)/(heat in). 
Here this is 1 — Qi/Qu . whence 



C p (T a -T b ) 

V C v (T c -T b ) 

Pg(Va ~ Vb) 

1 V h (P b -P a ) 

At first glance, there is something odd about this expression: we expect that the efficiency is 
dependent only on temperature, and independent of the properties of the working substance, yet 
C p and C v appear in the expression. The resolution to this is that we cannot pick the pressures 
and volumes independently: given the two state points at V bj the properties at V a are completely 
determined by the adiabat. 

10. Work and thermal equilibrium 

We have a Carnot engine, but the temperature of the two reservoirs is not a constant. So consider 
an infinitesmal amount of work done while the temperature can be regarded as constant: from first 
law dW = dQ ' - dQ" 


26 


Download free eBooks at bookboon.com 







An inverted textbook on thermodynamics 
Part II 


Cycles and the Second Law 


As heat is removed, the temperature of the heat bath drops (and vice versa), by dQ = CpdT, 
where Cp is the heat capacity. 

Now integrate over all temperatures as the heat baths approach equilibrium temperature Tf. 

r T f r T f 

dW = -C P ( / dT' + / dT") = C P (Ti + T 2 - 2 T f ) 

J T\ JT 2 


This still tells us only that some heat was converted into work. To calculate Tf we must apply the 
Second Law. 


For a Carnot engine, the entropy of the universe is conserved so that for each infinitesimal heat 
transfer, the increase in entropy of the low temperature reservoir equals the decrease in entropy of 
the high temperature reservoir. 

dQ' dQ" 


Taking advantage of the fact that the heat capacities of the two bodies are identical, and being 
careful with signs: 


r Tf dV f 7 

Jt i T' ,/t 2 


Tf dT" 

Jfff 


= In T = ^ TL 
T f T 2 


So Tp = VTiT 2 , and the total work done is = Cp(T\ + T 2 — 2\/TiT 2 ) 

This reflects the situation in any case where the amount of heat energy is finite: by doing work 
the engine returns the two bodies to equilibrium. Only by continually supplying energy can the 
reservoir temperatures be kept constant. 

If the cold reservoir is large, then we can assume the temperature T" = T 2 = T/, however the heat 
capacities are no longer equal. 

So the First Law gives: W = Cp(Xi — T 2 ) — f dQ" 



A cate-Lucent 


www.alcatel-lucent.com/careers 


What if 
you could 
build your 
future and 
create the 
future? 


One generation’s transformation is the next’s status quo. 
In the near future, people may soon think it’s strange that 
devices ever had to be “plugged in.” To obtain that status, there 

needs to be “The Shift”. 


27 

Download free eBooks at bookboon.com 
















An inverted textbook on thermodynamics 
Part II 


Cycles and the Second Law 


and the Second Law gives: 


From which we deduce: 



dV 

Yd 


/ 


dQ" 

T 2 


Q" = C P T 2 In F 
J-2 


W = C P T\ - C P T 2 ( 1 + In F) 

-L 2 

If we plot these expressions, we see that no work is extracted if T\ = T 2 . Otherwise, the second 
expression is always larger than the first. 

The second case is closer to the situation in a power station, where a finite source of heat is dumped 
to a cold reservoir. It illustrates the importance of having cold water, either direct from a river or 
via a cooling tower, for the cold reservoir. If the water in the “cold” reservoir is heated a bit, it 
can be used for district heating. 


5 Entropy 


1. Calculation of entropy change 


The heat flow must be sufficient to change the water temperature by 20K, i.e. Q = 5 x 20 x 4.19 = 
419 kJ. 


Entropy change is given by 


AS = 


/ 

J s: 


finish 

start 


dQ 

T 


For the surroundings, the temperature is a constant, 278K, so AS = Q/T = 419/278 = 1.507 kJK 1 . 

For the water, its not so simple because the temperature changes. We use the fact that entropy is 
a state variable, so the difference in entropy of the water is independent of the path between the 
start and finish of the process. 

So imagine a reversible slow heating process such that dQ = mc p dT , 


AS = 



mcpdT 

T 


i ( T M 
J-2 


5 x 4.19 ln(278/298) = -1.46/c JR- 1 


A very simple calculation, but a subtle assumption: we can calculate the entropy change in an 
irreversible process by considering an equivalent reversible process. We can do this because entropy 
is a state variable. Also, the heat Q flows out of the water which means that dQ should be negative: 
in the integral this appears because dT is negative in the cooling process. 

Note also that entropy is not conserved. The increase in energy of the surroundings is equal to the 
loss of energy from the hot water, but the increase in entropy of the surroundings is bigger than 
the reduction in entropy of the water. 

Finally, notice that in a process of heating water from 5°C to 25° C, the total entropy still goes 
up. The general function for the entropy change is f(pc) = x — ln(l + x), based on x = AT/T. 
This function is positive for any value of AT > —T. It becomes infinite at x = — 1, i.e. cooling to 
absolute zero would generate infinite entropy in the surroundings. 

2. Variation on the same theme 

The resistor is in a steady state, so it has constant temperature and internal energy. If we look 
at the resistor some time later, it is identical (still has electrons flowing through, so not strictly in 
equilibrium, but not changing in time. In fact, all its state variables are constant, so A S res i s t or = 0. 

Electrical work is done by the current, and this is converted into heat by the resistor. All this heat 
goes into warming the water: entropy is continually being produced, if we consider 1 second. 
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Since the temperature is constant, the entropy can be easily found with the normal integral 




/ 


dQ 

T 


2000 J 
300 K 


6.7 J/K 


Hence the entropy increase in the universe is. 

A*S total — A*S 'water H - AS resistor — 6.7 J/K/ SCC 


All the electrical energy is converted to heat at 300K. So an electrical heater has a 100% efficiency 
(77 = 1). Don’t be impressed: the electricity generation would have been less than 100% efficient, 
and a heat pump would have rj > 1. 

In reality, the temperature increases by a small amount, which would make the equation somewhat 
more complicated, but the entropy change must be bounded by all heat supplied at the input 
temperature Q/300 and all heat supplied at the output temperature Q/(300 + ST). 

3. Entropy and ideal gases 

Use central equation TdS = dU + PdV to write 



PdV 

T 


dU dV 

-t +r v- 


Now we use the fact that the energy of an ideal gas depends on its temperature only (the actual 
value of this energy varies between monatomic, diatomic gasses etc.) 

So for a constant T process, U is constant and so dU = 0 and 


AS 

- < 1 

(i) 

S 2 ~ Sy 

- 

(2) 


If we think of entropy as being the number of different ways of arranging the molecules, then it is 
intuitive that larger volume means larger entropy. 

Using the definition of Cp in terms of entropy Cp(T, V) = T (§fr) p , the pressure dependence of 
Cp(T,V) at constant T is given by 


dCp\ 
dP ) T 


= T 


= T 


d 2 S 


dPpdTpJ 
d 2 S \ 
dTpdPp) 


= 


\dT 2 


( 3 ) 

( 4 ) 

( 5 ) 

(6) 


where we used the fact that for a state variable, the second derivative is independent of the order 
of differentiation (commutative), and then used a Maxwell relation to eliminate S. So far this is 
just another example of a general relationship between materials properties. 
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For an Ideal Gas we use the equation of state. 


d 2 V 

Qf2 


d(R/P) 

dT 


Since R is always constant, and P is constant here. 


4. Entropy changes in the ideal gas 


A slightly more complicated version of the previous question: again we use central equation TdS = 
dU + PdV to write 



PdV 

T 


We would like to eliminate dU , which we can do by considering the first law dU = dQ — PdV and 
differentiating with respect to T. 


fdU\ _ f dQ\ _ p (dV\ 
\dT ) v \dT ) v \dT ) v 


Now by definition the first term is 


Cy = 


(dQ 

\0T 


v 


And trivially, (j^) v =0 since there is no change in volume at constant volume. Hence 

cv = (jVj ^dU(T,V) = CvdT + 0 


dS(T, V) 


C v dT PdV 

~T ^ T 
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So far this is all general for any material at constant volume. 

In the specific case of the ideal gas (defined by PV = nRT ) we can use the first equation to create 
two integrals with single variables: 


dS = C v 


dT 

~T 


T tiR 


dV_ 

~V 


AS 


r T2 A + BT 

JTi T 


dT + nR 



dV 

~v 



+ B(T 2 — T\) + nRln 



( 7 ) 

( 8 ) 


Now this should look worrying: AS appears to diverge as T\ goes to zero. This is still true for 
B = 0, the normal ideal gas. Could an ideal gas really have an infinite (negative) entropy? In fact 
quantum mechanics comes to the rescue - classically T=0 implies we know the precise velocity of 
the particle (zero) and V=0 implies we know its position exactly. This violates the Uncertainty 
principle, so the idea of an ideal gas breaks down at T=0. 

There is an expression called the Sackur-Tetrode equation which gives the “Entropy of an ideal gas”, 
unfortunately it involves approximations which means that it breaks down at low temperatures. 

5. Another look at the Carnot cycle 

On a TS diagram, the Carnot cycle of isotherms and adiabat/isentropes is just a rectangle. Define 
this as being between T 2 > T\ and S 2 > S 1 

Heat is absorbed only in the isotherms, an amount dQ = TdS , so assuming the path runs clockwise 
(Carnot engine), the total heat absorbed is: 



T 2 dS + 



T 1 dS=(T 2 -T 1 )(S 2 


Si) 


Which is the area inside the cycle. So the net Heat absorbed is the area inside the cycle in a TS 
diagram, while the net Work is the area inside the cycle in a PV diagram. And of course they are 
equal to one another. 

Mathematically, what we just did was a contour integral of dQ around the cycle. If we go the 
other way round the cycle we get negative heat in: a Carnot refrigerator. 

6. Proving that booze increases randomness 

The ice absorbs heat as it melts, but the water stays at 0°C. So the heat must be coming from the 
air. Consider first the ice, where everything is happening at 0°C and we are losing latent heat. 


ASi ce — AQi ce /Ti ce — 

J- ire 


300 x 334 
273 


1 


The entropy of the ice increases as it melts 

Now consider the water. It stays conveniently at 0°C so no net heat is absorbed by it. AQ water — 
—AQi ce + AQ a ir = 0 so J dS J dQ/T = 0 
Finally consider the air 

A S air = AQ air /T air = = ~ 3 °° * 334 = -34.2 J 

-L air Zb/O 

The air loses heat and its entropy is reduced. In fact, it should cool down a tiny bit. 

Overall the G and T increases the entropy of the universe by 2.5 JK~ X . 
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7. Entropy change, reversible and irreversible processes 

Divide the universe into the system (the lead, Pb), surroundings (heat bath) and everything else. 
Assume that the heat bath is thermally isolated from everything else, so the only entropy changes 
are in the lead and the bath. 

Dropping hot lead into water is a rapid, irreversible cooling process. However the change in entropy 
of the lead is the same in each case, because entropy is a state variable, and the lead begins and 
ends in an equilibrium state. So we can calculate it by considering a reversible route between the 
states, successive equilibria between a series of heat baths at infinitesimally different T. We know 
AQp& = CpAT so the entropy loss is 

AS Pb = J d Q r = J = C p ln[100/200] = -1000. In 2 = -693 J/K 

(a) We conventionally make the assumption that the bath is “large” and maintains its temperature 
at a constant, call it Ip, so its entropy is increased via an isothermal process: 

AS Mh = £2 = j iT = 1000 ( 200 - 100 ) = K 
Tb Tb J 100 ' 

The overall entropy change of the universe is then (1000-693) = 307J/K. 

(b) With two baths we use the same equations 


AC A Q , A Q 1000.(200 - 150) , 1000.(150 - 100) OOOT/ ^ 

ASbaths = ^ ^ = -Y50- + -100- = 833J/K 

The overall entropy change of the universe is then (833-693) = 140J/K. 

(c) In a reversible process the entropy of the universe is conserved. 

The question illustrates that the entropy of the universe depends on the process, even though 
the entropy of the system does not. The reversible case required a series of heat baths out of 
equilibrium with each other. An alternative would be to use the heat transfer to drive a Carnot 
engine, in which case some of the heat would be converted into work. This gives an alternative 
view of work: energy which does not contribute to the entropy. 
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6 Thermodynamic Potentials 


1. Heat Capacities 

For the heat capacity at constant volume, from the First Law dU = dQ + dW , differentiate with 
respect to T at constant volume: 


(dU_ 

V<9T 


v 



= Cy 


where work is — PdV , and dV , the change in volume, is obviously zero for a constant volume process 

Now using the Central Equation dU = TdS — PdV. Once again dV is zero for a constant volume 
process, so: 


Cy 



= T 



v 


For the heat capacity at constant pressure, directly differentiating the First Law would give two 
terms. 


Cp = 





p 


More convenient is to first write the First Law in terms of enthalpy (H=U+PV; dH = dQ+VdP). 
Now for a constant pressure process, dP=0. 




Finally, using the Central Equation, dU = TdS — PdV => dH = TdS + VdP. So 


( dH 

VST 


p 



= C P 


It is not essential to introduce enthalpy, but it makes the analysis of constant pressure processes 
much more straightforward. 

2. Helmholtz function and pressure 

Helmholtz free energy is F = U — TS, so using the central equation: 


dF = dU - TdS - SdT = -PdV - SdT 


Expanding F = F(V,T) we can write 


so we can equate 


dF = 



dV + 

T 



dT 



S = - 



v 


So far this is general. Now for this particular gas: 


P 



— RT\tl(v — b) + /o(T) 


a RT 
v 2 v — b 
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Where fo(T) is some constant of integration, which does not appear in the equation of state. 

Notice that we can tackle these problems using either total free energy, entropy, volume etc. (F, 
S, V), or specific quantities (f, s, v). Pressure and Temperature are the same in each case. 

3. Using Maxwell relations 

By definition F = U — TS, so: dF = dU — TdS — SdT , which combined with the central equation 
dU = —PdV + TdS , gives 


F = -PdV - SdT 


So we can equate 



We use the fact that for a state variable (F) the order of differentiation does not matter. Considering 
second differentials gives: 



Which leads to a Maxwell relation: 



To get the volume dependence of c v =Tjj^ v we find: 



Now we are given an expression relating the heat capacities 



and by definition 



p 



Now from the central equation, and application of the Maxwell relation from above, gives 



so that, using the cyclical relation 



(ii) 


( 10 ) 


( 9 ) 


VTP 2 

K 


( 12 ) 

( 13 ) 


34 


Download free eBooks at bookboon.com 










An inverted textbook on thermodynamics 
Part II 


Thermodynamic Potentials 


Now look at ratios. By definition, the compressibilities, Pp = — y (§^) T ; p s = — y 


cp 

c v 



(14) 


we have (ff ) p ® s (§f ) T = -1 and (§f ) v (§£), (f ) T = -1, so that 


(dT\ (dV\ (dV\ 

C P _ \dV) S V dS ) T _ \dPJT _ K P_ 

^ “ (§?U§§)t“ Ws'** 


(15) 


Using the results Cp — Cy = and 

We can write 


Cy Pp 

PT ~ PS = PT{ 1 -) = « T (1 - 7 —) = 77 “ 

ft T Op Op 


VT /? 2 


Pp 


VT(3 2 

Cp 


This is all general. For the specific case of one mole of ideal gas: PV = RT, P ( 'jyj-) p = R\ 
(3 = 1/T; K = 1/P; P (fp) T = —V/P. So that 


c p -c v = VT(l/T 2 )P = PV/T = R 
where we use the cyclic relation (§f ) v (f^)p (§p) T = —1 

k t -k s = VT(1/T 2 )/^R = 2V/5RT = 2P/5 
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Notice that for a diatomic gas the result would be different: kt — ks = 2P/7, and for a giant 
molecule (or solid) the two compressibilities become essentially the same. Tables of compressibil¬ 
ities and other elastic constants for condensed phases seldom distinguish between isothermal or 
adiabatic. 


4. A block of metal. 


“Adiabatic” means no heat flow dQ = 0, which in the reversible case implies dS = 0. Choosing 
variables S(T,P ) lets us write 

ds= (§) P dT+ (w) T dp= ° 

Using the Maxwell relation (§p) T = — (fyr)p = —V/3 and the definition Cp = = T (ff ) P 

it follows that 

S vfdp= ! 


= I ^fdT 


Solving the integrals, assuming volume is constant 


whence 


V(3[P 2 - P,} =c P In ^ 

J-l 


1 T2 V P \P PI 

= - [P 2 -P 1 ] 

±1 Cp 


This question is of a classic type beloved of thermodynamicists: “what happens when you change x 
and y, holding 2 constant”. The technique is to expand one variable in the other two, eliminate the 
term in dz, then use the Maxwell relations and/or triple product along with definition of materials 
properties to get to measurable quantities. There are several ways to get the answer, so lets do it 
again: 

Expand T(5, P): 

'dT N . 

Jp). 


dT=l "is) P dS+l ^ dF 


For an adiabatic process, dS = 0, so: 


Triple product, then Maxwell relation 


dT 




dS 


dS 


dP 


dT=- — — dP = —- — dP 


dT 


dS } p \dT 


dV 


Identify materials properties: 


dT = —Vf3dP 

Cp 


Rearrange to obtain the same integrals as before 

/ 

and the correct answer follows once again! 


T /3 
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5. From Gibbs function to equation of state 

We are given the specific Gibbs Free Energy 


g = RT\n P + A + BP + -CP 2 + -DP 3 

z o 


Expanding g = g{T , P ) we can write 



T 


dP 


and by definition (g = u + Pv — Ts ) we also have dg = vdP — sdT , so that we can identify 


V 


T 

DT 

— + B + CP + DP 2 



which is the equation of state. The constant A simply sets the zero of energy, which can be chosen 


arbitrarily. Often one can measure the equation of state experimentally, or calculate the Gibbs 
function theoretically. It is useful to know that, given one, it is possible to obtain the other. 

6. A harmonic material The equation of state is given by P = A(y — b) + CT 

We are given the form of the equation of state, but not the parameters. To determine them, we 
must build a set of simultaneous equations, one from each of the given properties. 

Taking the value of v and T at p=0 we find Bulk modulus K T = 10 10 = —Vo (§y) T = — 10 -3 A; 
So A = —10 13 . 

Expansivity 0 = MItOp = £ ( B)t (§?)„ = 10_5 ' 

From which (3K T = 10 5 = (§f ) v = C 

Now use the pressure: 0 = A (lO -3 — i;q) + 300C; 

Rearranging = 300 C/A + 10 -3 = 1.003 x 10 -3 . 

Inspecting the equation of state, there is a non-infinite equilibrium volume at P = 0, even for 
high-T. Something must be holding the material together - this is a condensed phase. 

The problem is oversimplified in that the measured properties should be pressure and temperature 
dependent. However for small changes in T and P the linear equation of state should be reasonable. 

7. Deriving the ideal gas equation from experimental laws. 

We start with the empirical observations known as Joule’s Law U = U(T) and Boyle’s Law PV = 
/i(T) where /i is some unknown function. 

Joule’s Law implies that (^) T = 0 for any quantity A, so we start with the central equation 
dU = TdS — PdV , and differentiate by volume at constant T: 



rearranging, and using a Maxwell relation gives a relationship for constant volume process. 



dP _ dT 
~P ~ ~T 
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Integration gives 


lnP = lnT + f 2 {V) 

P/T = f 3 (V) 

So P/fs(V ) = T, and Boyle’s Law gives PV = /i(T). The only way to make these compatible is 
if fz(V) = ai/V and fi(T) = a 2 T, with undetermined constants cq. So that 

PV = aia 2 T 

and we can identify the constants with the gas constant: aia 2 = nR. 

Once again, the equation of state can be deduced from experimental observations. 
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7 Expansion Processes 

1. Free expansion of the van der Waals gas 

Free expansion conserves ?7, because there is no heat input and the expansion is done into a vacuum 
(i.e. zero pressure), so no work is done. Starting with the Central Equation dU = TdS — PdV , we 
get an expression for a volume derivative at constant T: 



We are trying to obtain an expression for and we can get this quantity into the expression 

by using a cyclic relation (f^) T = - (§^) y 





-P 


Using the definition of Cy = (j^j = (§^)y? we get 


-a 




-p 


Entropy does not appear in the van der Waals equation, so we must eliminate S. This can be done 
using a Maxwell relation: (§y) T = (§^)y> so substituting and rearranging 



- pyc v 


So far this is completely general. Now consider the special case of the van der Waals gas: 


RT a ( dP \ _ RT 

V^b~V 2 ’ \df) v ~V^b 


(dr 


u 


~[T 


R 

V-b 


RT + a ]/r 

V^b + y 2 \/Cv 


—a 

CyV 2 


Note that for an ideal gas a = 0, so the temperature doesn’t change in a free expansion, however 
for a non-ideal gas the important term is a, which represents the interactions between atoms. The 
sign of the cooling/heating effect depends on whether the gas atoms attract/repel. This is because 
the internal energy is split between bonding (potential energy) and kinetic. If the potential energy 
is increased (longer bonds), while V is constant, then the kinetic energy must be reduced. 

The van der Waals equation of state works well for Helium at room temperature, with a > 0 and, 
so we see that helium cools when it undergoes a free expansion at room temperature 

2. Joule-Kelvin coefficients 

(a) This is another piece of manipulation of partial derivatives. We will need to absorb H into 
the heat capacity, so, start by taking triple product 

fdT\ fdT\ fdH\ 

\dPj H ~ {dH) p {dPj T 
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We also know from dH = TdS + VdP, that 




and 


dH 

~dP 


= T | — 
1 dP 


dV 

V = —T — 

1 dr 


V 


where we used a Maxwell relation to eliminate the entropy. Now, using the definition of the 
coefficient of thermal expansion we have 




dP 


H 


Cp 


(b) For an ideal gas, PV = nRT we have 


dV 

a ~V [dT 


1 nR _ 1 
_ T 


So that (aT — 1) = 0, and hence /ijk = 0, so an ideal gas doesn’t change temperature in an 
isenthalpic expansion. 

(c) For the van der Waals gas, it is a bit fiddlier. To get the thermal expansion, it’s easiest to 
write the equation of state in terms of T, using specific heats and volume: 

RT = [P + ^}(v - b) 

V z 


and then find: 


R 



(P + ■%) - 2 a(v - b)/v 3 
v z 


RT _ 2 a(v - b ) 
(v — b ) v 3 


Where we eliminated P so that the expression refers only to two variables (T, v) which can be 
taken as independent, so that: 

v f RT / v 

^ JK cp \RT/(y — b) — 2 a(v — b)/v 3 

Remember that the specific heat cp is not uniquely determined by the van der Waals equation. 
It depends on whether the gas is a monatomic, diatomic, etc. So we can do no better than 
leave it as cp, and remeber that it must be positive. 

(d) The inversion curve is where hjk = {§p) H = 0? i- e - when the effect of a Joule-Kelvin 
expansion changes from cooling to heating. Therefore, from setting the bracket to zero in the 
expression for /jljk, the inversion curve is specified by: 


—— = RTin/(vin ~b) - 2a(viN ~ b)/vj N 

VJN 


which after some tidying up gives an equation for the inversion line 


RTin = 


2 cl(vin — b ) 2 


IN 
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Sketching the curve shows that the highest value for Tjjy = 2a/b comes at the high-?; asymp¬ 
tote. Above this temperature, Joule-Kelvin expansion of an ideal gas always results in heating. 
To obtain this formally, we need to differentiate to find the turning point: 


dTny 

dviN 


2 a 


(vin ~ b) = 0 


Mathematically, this implies a turning point when vin = b. Plugging this back into the 
equation for the inversion temperature gives a value for zero. This is a minimum of the 
function, although v < b is unphysical as it implies that the volume available for the substance 
is less than the volume of the component atoms. 


Now for some mathematical legerdemain. Rearranging the inversion line equation and square 
rooting gives 


RbTiN 
2 a 


(vin - b) 

VIN 


from which 


2 ab 

VlN = 2 o^rFTin 


which we can substitute back into the van der Waals equation to get an expression for the 
“inversion curve” in PT space: 


P 


8a 

2762 


9 — 12( 


8aT 
27 Rb 


\/3) 2 


This gives a maximum inversion pressure at T = Above this pressure, the Joule-Kelvin 

process always results in heating. 
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3. Cooling in an adiabatic expansion 


Figure shows curves P=l/V, P=l.l/V, PV 7 / 5 
representing two isotherms and an adiabat 
(with nR=l). It can readily be seen that be¬ 
tween the isotherms the adiabat has the larger 
slope. The temperature at any point along the 
abiabat is that of the isotherm passing through 
that point. Hence expanding the ideal gas adia- 
batically takes it from one isotherm to another 
corresponding to lower temperature. 



For the adiabatic expansion process we are interested in the change in temperature when pressure 
is reduced at constant entropy (no heat exchanged) (f^)s'- 

To get this into a managable form, we eliminate S using the triple product, then a Maxwell relation, 
then the definitions of Cp and (3: 

f dT \ _ f dS \ f dT \ -f dV \ f dT \ -f dV \ T _ TVf3 
\dP ) s ~ ~ \dP ) T \dS ) P ~ \dT ) P \dS ) P ~ \dT ) P C^ ~ 

So the adiabatic expansion always produces cooling, and the cooling is always larger than for the 
isenthalpic process: (§p) H = TV q~ V • 

For completeness, the coefficient for the Joule process: 

_ (dT\ _ (dT \ ( dU \ 1 

~ WJh - ~ \du) v Vdv) T ~ c~ v 

has different dimensions, so cannot be directly compared. 

4. The Joule-Kelvin process involves a continuous flow of gas which can be achieved in is cyclic 
process. By contrast, the Joule process involves a one-off expansion into a vacuum and cannot 
easily be used for continuous cooling. Neither process is reversible, so they cannot be used in a 
Carnot refrigerator. The adiabatic expansion process lowers the temperature of the working fluid 
by more than the Joule-Kelvin process, but since it doesn’t extract heat from the environment it 
is not useful in refrigeration. 

5. Expanding through a Phase Transition 

To maximise /ljk we want a very large value of In a boiling phase transition, the volume 

changes discontinuously at the transition: a is infinite, as is Cp. So the equation implies integrating 
over a divergent function. 

6. Gibbs Helmholtz Equation 

From the definition dG = Vdp — SdT , 

Now consider the quantity: 

dG/T\ _ 1 fdG\ G _ (G + TS) _ H 
dT ) p ~ T [dT) p ~ - Y 2 “ ~T 2 

Which in integral form can be written 

I d{G/T) = — J Y T 


we have 


fm - 


\dT J 


P 


= -S. 


dS- 

dV 


-P 


T[3K t 

C v 
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doing the integral gives... 


and rearranging 


G 2 _ G i 




g 2 


GiT 2 

Ti 


T 2 


f 

JT ! 


t 2 h 

—rdT 
2^2 


We now apply this equation to two different phases of the same material. Let state 1 be atmospheric 
pressure and 300K, state 2 atmospheric pressure and 400K. We are told that 300K is the melting 
point, so here the gibbs free energies are equal: G\ = Gf = Gi, and the free energy difference 
G l - G s is zero. 

At 400K: 

A G 2 = G\ - G'l = -T 2 Hl ~ 2 Hs dT 
Calling the volumes V\ and V s we have 



§■dT = [Ui + PVi]( 1/Ti - 1/T 2 ) + A, ln(T 2 /Ti) 


^dT = [U s + PV S ](1/Ti - 1/T 2 ) + ^ ln(T 2 /Ti) 


So that at T 2 

AG = -(T 2 /Ti - l)[tf z - + P(y z - V 3 )] - (A z - A s )T 2 In(T 2 /Ti) 

(sanity check: AG = 0 for T 2 = Ti) 

The first term refers to the enthalpy difference between the phases. We can expect the constant 
[Ui — U s + P(V/ — T4)] w iH t> e positive, since the solid has lower enthalpy that the liquid. We can 
also expect (Ai — As) to be positive: these quantities are heat capacities, and we expect the liquid 
to be higher. Taken together, AG is negative for T 2 > T\ and positive for T\ > T 2 - the stable 
phase of the material is liquid above the melting point and solid below. 

As T —> 0 the free energy difference tends to the enthalpy difference, as it should. However, with 
the given expression the heat capacities remain finite, in violation of the Third Law. 

Notice that we do not have enough information to calculate the actual values of G, just the free 
energy difference 

In analytic thermodynamics, the Gibbs-Helmholtz equation tends to look like just rearranging 
symbols. In fact, it is incredibly useful in Monte Carlo and Molecular Dynamics calculations which 
are the mainstay of computational physics. Given a set of atomic positions and velocities, it is easy 
to calculate the enthalpy directly from the Hamiltonian. However, there is no way to obtain the 
Gibbs free energy, or anything involving entropy. To understand why, remember that the entropy 
represents all the possible ways the atoms could be arranged to give the same macroscopic phase. 
Knowing what a typical state looks like tells us nothing about how many there are. The Gibbs- 
Helmholtz equation gives us a way to determine the changes in entropy (number of states), simply 
by measuring changes in enthalpy (energy of states). 

Notice that it does not allow us to determine the value of G, only values for changes, AG. 
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7. Critical point 

The van der Waals gas has P = RT/(v — b) + a/v 2 . 

Considering turning points with = 0 we find that there are three cases for an isotherm 
T > T c : no turning point on the isotherm with = 0. 

T < T c : one maximum point on the isotherm with ^y = 0 and < 0. 

T = T c \ one inflection point on the isotherm with ^y = 0 and = 0. 

there is also an unphysical region with v « b where either the pressure or temperature must be 
negative. This region has an unphysical negative bulk modulus. 

With T < T c , for given external pressure and temperature, there are two possible volumes. What 
actually happens here is that if the external condition fixes the volume, we get a coexistence 
between the two phases, in the right proportions to fill the volume available. Alternately, if the 
external pressure and temperature are fixed, any infinitesimal difference from the values on the 
phase boundary causes the material to transition int the relevant phase. 

To find T c we require both of the first and second differential of pressure to be zero. For a van der 
Waals gas, 


P = 



nRT an 2 
V -nb ~ W 

nRT ( 2an 2 

(V - n6) 2 + 

2 nRT 6an 2 

(V - nb) 3 ~ 


= 0 

= 0 


Eliminate T from the above: 


nRT 3 an 2 3 an 2 3 nRT 

(V - nb) 3 = ~W~ = 2VW = 2V (V — nb) 2 

From the first and fourth expression 2V = 3(V — nb); V c = 3 nb. 

Remembering the idea that b is the volume of the atom, this is when “The atoms take up one third 
of the volume”. 

By substitution, we can find T c 

2 an 2 (V c — nb) 2 2 an 2 4n 2 b 2 8a 

c V c 3 nr 27 b 3 n 3 nR 27 bR 

and finally 

nR(8a/27bR) an 2 

c 3nb — nb (3 nb) 2 

4 a a 

27h 2 ~ 9b 2 
a 

27tf 
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Sometimes, it is convenient to treat the van der Waals gas using dimensionless units such that the 
triple point lies at T c = 1, P c = 1, V c = 1. This saves us having to carry round factors of i?, a and 
6, at the cost of not being able to check dimensions of our equations. 

Supercritical fluids The planetary atmospheres of Venus (mainly CO 2 ), Jupiter and Saturn (mainly 
H 2 ) are supercritical fluids. They combine the liquidlike ability for high solubility with gaslike high 
diffusivity, so can be used e.g. in decaffeination. They are also immune to microcondensation and 
bubbling. The combination of high heat capacity with high diffusivity, makes supercritical water 
an ideal coolant in power stations, where heat must be transferred efficiently to the generator from 
the reactor/“boiler” (obviously “boiler” isn’t quite the right word here!) 
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8 Thermodynamics in other systems 


1. A rubber band 

a) Start with Central Equation dU = TdS — PdV from which, differentiating by V. 

=t(™) -p 

dV) T \dV) T 

Use Maxwell relation (§y) T = (§§) v to get the required answer 


dU 

W 


= T 


dP\ 
dT ) 


-P 


v 


b) For a wire, substitute P (pressure) with -F (force) and V (volume) for L (length). Note that 
increased pressure makes the volume smaller, while increased force makes the wire longer. 


&U\ = ^ T fdF\ 

dLJ T \dTJ 


) 

T V / L 

You might think the second term looks odd, normally force is just the derivative of the energy? 
What’s happening here is the boundary conditions, we’re not simply stretching the band, we’re 
also allowing it to absorb/emit heat. The normal spring definition is recovered if we consider an 
adiabatic process on dU = TdS + FdL : 

'dir 

~8L 


= F 


Remember that F is the force exerted on the band, not by it. 
If we have an isothermal process, then rearranging the above: 

F= (m +r( a -L 


dL J T \dT 

^ -t(—\ ( 9 ( u ~ ts ) 


dL 




dL 


and using a Maxwell relation 

F= 'dL) v _ , 

m J J rj~i y I J J rj~i y # J J rj~\ 

We recover the result that the required force for an isothermal extension is the derivative of 
Helmholtz free energy , U-TS. 

c) The equivalent of an equation of state for a wire is F = aT — (^f) 2 
From which we can calculate 


dU 

dL 


= -aT 


L 

Lo 


aT 


L 

Lo 


= 0 


i.e ?7(L, T) = U{T ): internal energy depends on temperature only. This works for this carefully 
chosen equation of state which has energy proportional to temperature. 

This does not mean that we can stretch the wire at no cost. If we stretch the wire to length L 
isothermally, the relevant thermodynamic potential is the analogy of Helmholtz free energy. 
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d) Work done, again by analogy 


W= [ aW = [ 2 FdL 


Jl x 



Using the values L 0 =lm, Li=lm, L 2 =2m, a=1.3x!0 2 NK x , T=300K. 


W = (300).(0.013)(2 - 0.5) + (300).(0.013)(0.5 - 1) = 3.9J 


Since V = U(T ) only, and T is constant, U is constant, so from the first law dU = dQ + dW = 0 
the work done is equal to minus the heat input. 

So heat input is -3.9J, alternately 3.9J of heat is given off. 

2. Adiabatic rubber band 

Now consider the adiabatic case, where we don’t give the wire time to warm up after stretching it. 
Adiabatic and reversible means dS = 0, so write an expression for the entropy. 



The relevant specific heat, at constant length, is cl = T (ffr) L , and we have a Maxwell relation: 

(ds\ --(9F) 

V dL ) T ~ V dT ) L 

So we have a relation 



We now integrate this equation 




Plugging in the numbers 


^ ln ^ = 0 - 013 -[ L5 -°- 5 ] 


T 2 = 300exp(0.013/1.2) = 303.3A 


Change in temperature is a 3.3K increase. 
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3. Entropy of diamond 

Given that cp = 124(T/1860) 3 . 

To find entropy change, consider a reversible process involving small changes at constant volume: 
(IS = c v dT/T 


AS 



124 r300 3 4 3 


I860 3 [ 3 3 

0.173 kJ/K/kg 


For graphite, similarly 
!&*[*£-£] = 0-237 kJ/K/kg 

Thus graphite has higher entropy and is more stable at high temperatures. Since it is also the less 
denser phase, it means that As and Av for the diamond-graphite transition have opposite signs, 
so the Clausius-Claperon slope is negative. 

As T —> 0 Cp —> 0. Cp involves changes in entropy, so the Third Law is satisfied. There is no 
divergence in the integral at T=0. 



Brain power 


By 2020, wind could provide one-tenth of our planet’s 
electricity needs. Already today, SKF’s innovative know¬ 
how is crucial to running a large proportion of the 
world’s wind turbines. 

Up to 25 7o of the generating costs relate to mainte¬ 
nance. These can be reduced dramatically thanks to our 
(^sterns for on-line condition monitoring and automatic 
lul|kation. We help make it more economical to create 
cleanSkdneaper energy out of thin air. 

By sh?fe|ig our experience, expertise, and creativity, 
industries can boost performance beyond expectations. 

Therefore we need the best employees who can 
kneet this challenge! 


Power of Knowledge Engineering 


Plug into The Power of Knowledge Engineering. 
Visit us at www.skf.com/knowledge 
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4. Planck’s Law 


We are given that P = U/3V. 

From the central equation dU = TdS — PdV , 



and using one of the Maxwell’s relations, 



Energy density is defined as u = U/V , and or the photon gas P = U = uV and u = u(T ), so 
we get 


u = 


It— 

3 dT 


1 

3 U 


which implies that 4^ = ^ Which integrates to 4 In T +In (^fr ) = hi u with appropriately defined 
for some constant of integration. 



The Planck distribution law, for energy density of cavity radiation as a function of wavelength A, 
and temperature T is: 


u\(X,T) oc 


1 

V 


1 




To get the total energy density, we need to integrate this over all wavelengths: 




ghc/ XksT 


dX 


If we make a substitution: 
x = hc/XksT then this becomes 


u = 27 rc 2 h 




—he \ 
k B Tx 2 ) 


dx 


Rearranging to make clear the dedimensionalisation the integral 


u = 27 rc 2 h 



x 


3 


e x -1 


dx 


The integral is clearly bounded, so it is simply a number. You can solve it if you enjoy that sort of 
thing (its about 27r), but the physically important result, the T 4 dependence, is thereby shown. 

Had we used the Boltzmann energy distribution, the integral would have been 




Which also gives a Stefan’s law T 4 result. It’s only if we use classical equipartition that the integral 
becomes infinite 


u = 


J u\(X,T)d\ oc J A A d\ 


oc 


the so-called “ultraviolet catastophe” since the integral blows up for short wavelengths. 
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5. Photon Gas 


(a) To obtain U(V,T) we will need to obtain a relationship between the internal energy, pressure, 
temperature and volume, then use the equation of state to eliminate the pressure. Start with 
the central equation 

dU = TdS - PdV 

Differentiate wrt V at constant T 



and use a Maxwell relation to eliminate S 



This is a general thermodynamic relation. We now apply it to the case of cavity radiation. 
Since specific internal energy u (energy density: energy per unit volume) depends only on 
T, the total internal energy depends on how big the cavity is, i.e. it is proportional to V: 
U = Vu(T) 



= u 


Using specific energy and volume, and P = ^ we find: 


T ( du\ _ 4 u 
3 [df) v -~3 


Rearranging and integrating: 




Integrating, and introducing a constant of integration In k 


In u = 4 In T -f In k 


Implies u = kT A or U = kVT A . 

(b) For a photon gas to have sufficient energy to move a piston: 

P = 10 5 Pa = u/3 = kT 4 /3 => T = ( 3P/k) 1/ 4 = 1.4 x lO 5 ^ 

This is hotter than the sun, but considerably less than in a nuclear bomb. The blast from a 
nuclear weapon is primarily a shock wave of air compressed by black body radiation. Since 
the pressure rises with the fourth power of the temperature, it quickly becomes enormous. 

(c) Since the energy depends only on T, the heat capacity is (^) v = 4 kVT 3 , or specific heat 
capacity 4 kT 3 . The Third Law requires heat capacities to go to zero as T —» 0, which is 
obviously true for a T 3 dependence. 

The specific entropy can be obtained considering a constant volume heating process from 
T=0. Using the Central Equation 

du = Tds — Pdv 

ds = ^'PdT = [ AkT 2 dT 
T dT J 

_ 4 kT 3 _ 4 u 

S ~ 3 “ 3T 
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(d) The Gibbs Free energy is defined as G = U + PV — TS , whence 




At equilibrium, only processes for which dG = 0 can occur. Remarkably, creation and absorp¬ 
tion of photons are both thermodynamically allowed because neither cause any change to the 
free energy. 

(e) We have derived G = 0 from the equation of state. However, from observations of the 
spontateous creation and absorption of photons, we could have used this as a starting point. 
Start with the definition 

dG = VdP - SdT 


since G is a constant (zero!) for the photon gas, dG must also be zero, so that: 


VdP = SdT 


Rearranging gives: dP = ydT Comparing this to the general expression for pressure as a 
function of T and some other variable X 


dP = 



dT + 



dX 


We see that (f^) T = 0, i.e. pressure is independent of any variable other than T. 
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6. Temperature of the sun. The total radiation from the sun is (47rR|)crT|, with a similar equation 
for earth (AttR^ctT^. 

The fraction of the solar radiation hitting the earth is given by the area of the earth divided by a 
sphere the size of the earth’s orbit (D): 

l*Rl) 

4ttD 2 

At equilibrium, the radiation hitting earth is the same as that emitted: 

(47I\R| )crT| = 

whence 


T s 


T e 


U D 2 




1/4 


2874 


/2.(1.5 x 10 11 ) 


(6.96 x 10 8 ) 


5959K 


This assumes that the earth’s radiation comes from the surface. In fact, the Greenhouse Effect 
means that much of the radiation comes from the (cooler and larger) upper atmosphere. It also 
neglects the contribtion of radioactive decay in the Earth’s core to warming the planet. 

7. Black hole entropy 

Start by writing the central equation as an expression for entropy 



PdV 

T 


Remember that we can write any state variable as a function of any other two: this time we use 

S(U,V) 

giving us (§§) v = So we need to calculate the relationship between S and U. 

We are given that S = ksAc 2 */AGTi, and we know that A = 47rr 2 ; r = 2 Gm/c 2 and U = me 2 , from 
which we find that S = so 


dS _ 1 _ Stt ksG 
dU = T = ftc 5 


A quick rearrangement gives 


Tic 3 

87 rksGM 


= 1.2 x io n A: 


Which is pretty hot, and so will evaporate very quickly. By contrast, a solar sized hole with 
m = 10 30 has T = 6 x l(T 8 iG 

In previous examples, we have been considering entropy in terms of atoms, however the Second 
Law applies to all processes, with increasing entropy determining what can and cannot happen. 
There are a number of difficulties with entropy in astrophysics, arising from the fact that a self- 
gravitating system is denser in the middle, and therefore doesn’t meet the equilibrium criterion of 
homogeneity. 
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9 Phase transitions 

1. Phases and Clausius-Clapeyron 

a) Sketch shows solid, liquid, gas, phase boundaries triple and critical points. 

b) In phases 1 and 2 at T,P coexistence, specific Gibbs functions are equal g\ = # 2 - 

dg = vdP — sdT 

at T, P. 


( dg\ 


( dg\ 


s \&r ) p ; v \dP) T ' 

At another point infinitesimally further along the coexistence line (T+dT, P+dP), still have g\ = (J 2 
which we can get by Taylor expansion from (T,P): 


gi (T + dT,P + dP) 

^+{%) T dP+ {w) P iT 

— s\dT + v\ dP 
(s 2 - Si)dT 



g 2 (T + dT,P + dP) 

^ T - p) + {w) T dP+ {w)/ T 

—S2dT + V 2 dP 
(V2 - Vi)dP 
(S2 ~ Si) _ l 
(v 2 -v i) T(v 2 -V i) 


Which is the Clausius Clapeyron equation. 

c) For the pressure cooker example, we plug in numbers to get 



p.t. 


A P 
AT 


= 3600 Pa/K 


therefore, for AT = 30K, A P = 1.09 x 10 5 Pa, so since we start at atmospheric pressure, the total 
pressure is 

P = 1.09 x 10 5 + 1.0 x 10 5 Pa = 2.09 x 10 5 , 

This is typical for a pressure cooker. “Cooking” is an irreversible process in which proteins are 
denatured - normally this means taking them from a non-equilibrium state to a lower free-energy 
state. Although we normally cook with temperature, it is also possible to use pressure alone to 
denature proteins. 

2. Phase mixtures 

Write the total mass of liquid and vapour as m = mi + m v . The mass is conserved in the vapouri¬ 
sation process. 

Total volume = V = Vi + V v == mivi + m v v v where v>i are the specific volumes (inverse density). 
Also define the specific volume of the whole system 


_ V _ vimi + v v m v 
m mi+ m v 


Rearranging this gives the “Lever rule” 


mi{v - vi) = m v {v v - v) 
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The name is by analogy with finding the balance point (pivot) of a bar with two weights hanging 
from it, m^m v being the weights, and v — vi, v v — v the distances from the pivot. 


The gibbs free energy of two phases in equilibrium is the same, and since neither T nor P change 
in the coexistence region, it remains constant. 


Compressibility is defined by 


P = 


1 dV 
V~dP 


and so, since the slope of the graph is zero, substances are infinitely compressible in the coexistence 
region. 


3. Solid-solid phase transitions 

The question asks about the phase transition line for a first order transition, so we will apply the 
Clausius-Clapeyron equation. Assuming dP/dT = AP/AT: 


AT = A P/ (— — -- ) = 291 x 99 x 10 5 x (- — t- ) /18500 = -5.8 K 

' \T(v 2 - Vl) J V 7333 5750 ) 

Like water, the high temperature, metallic white tin phase is denser than the low temperature, in¬ 
sulating grey tin. Consequently high pressure favours denser the white tin phase, and the transition 
temperature drops. 

4. Triple Point 

a) At the triple point, the vapourisation and sublimation curves meet, so we must have InP = 
0.04 - 6/T and InP = 0.03 - 4/T 



qaiteye 

Challenge the way we run 


EXPERIENCE THE POWER OF 
FULL ENGAGEMENT... 


RUN FASTER. — p 

RUN LONGER.. 

RUN EASIER... > 
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Whence Tt p = 200 K and by substitution: 

P tp = exp(0.03 - 4/200) = exp(O.Ol) = 1.01005 

b) To get the latent heats, we use the Clausius-Clapeyron equation: 

/ dP pb \ _ l _ l 
\dT pb ) T(v2~vi) Tv 2 

where the suffix pb reminds us that this is the derivative of the phase boundary line. 

We now assume that the gas phase V 2 can be taken to be an ideal gas. This lets us eliminate 

V\ « V2- 


Integrating this: 


dP _ l _ IP 
~dT ~ Tv2 ~ RT* 



InP = 


-l 

RT 


const 


Which can be compared directly to the given equations for the sublimation and vapourisation 
curves to give l = 6R and l = 4 R respectively. 

c) Consider reversible changes in a loop around the critical point. The entopy change around the 
loop AS = 0, since entropy is a state variable. 




For a loop close to the triple point, all temperatures are the same, so noting that e.g. Isv = — Isv, 
we get 


Lsl = Isv — Ilv = 2 R 


5. Impossible Phase diagram 

(a) The stable phase minimises the Gibbs free energy, dG = -SdT+VdP, so at constant T, increased 
P will favour smaller volume. Hence the density increases in order pi < /?2 < P 3 < Pi- The 
compressibility of phase 1 must be large, since at low pressure is has the lowest density, and at 
high pressure the highest. 

(b) For first order transitions, the denisties are different. At the triple point, we have pi < p 2 < 
P 3 < pi ~ phase 1 cannot have both highest and lowest density at the same (T,P). Thus the diagram 
is impossible. 

(c) Z = T — mP , which is a state variable. Consequently every point along the iso-Z line has a 
different value of its conjugate variable, Y which will be some combination of s and v. Precisely 
what isn’t important, all that matters is that this conjugate variable increases monotonically along 
the iso-Z line and changes discontinuously at the phase boundary. 

(d) At the triple point, since a first order transition implies the diagram implies that, Y\ < Y 2 < 
T 3 < Y <1 which is clearly impossible. 

This is known as the 180° rule, and is an effective shortcut for detecting incorrect claims regarding 
phase diagrams. It does not prohibit a maximum in phase boundary (e.g. the left hand figure if 
phases 2 and 3 were actually the same phase). At the point of maximum temperature, the densities 
of the two phases must be the same, at high pressure, phase 1 has higher density, while at low 
pressure phase 2/3 is denser. 
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6. Pressure effect on melting ice 

Another Clausius-Clapeyron question, and again assume the phase boundary to be a straight line, 
so that dPjdT = AP/AT 


A P 
AT 


l 


T{v 2 - Vi) 


333700 


273 x ( 10Q0 gi6 ) 


-1.3 x 10 7 Pa/K 


Now the applications 

(a) Pressure at bottom of Glacier = pgh = 916 x 9.8 x 2000 = 1.79 x 10 7 Pa 
Change in melting point: 1.79 x 10 7 / — 1.3 x 10 7 = —1.3 8K. 

So the melting point is —1.38° 

at the given temperature, —1.5° the ice is frozen. 

note we assumed it was frozen by using the density of ice (916). that assumption is now validated. 
There a small curiosity that if we assumed the density of water, the pressure would have been 
higher and the conclusion would be that the bottom was liquid. 

(b) Using the given value of —22°C, and assuming atmospheric pressure is negligible 
AP = (-1.3 x 10 7 ) x (-22) = 2.8 x 10 8 Pa 

which is about 2860 atm., validating our assumption. Above this pressure the crystal structure of 
ice changes to an arrangement denser than water. (The “Ice Ih - ice III” transition) 

(c) Skate has contact area of (0.3m) x (0.001m) 
assuming a 70kg skater, he exerts a force mg = 70x9.8N 
Increase in pressure is AP = ^ = Q 3 Q x X q 9 qqi = 2.3 x 10 6 Pa 

So A T m = = —0.18-K” 

This is too small to melt the ice in most cases. Even if we sharpen the blade to 10 _4 m, the drop is 
only 1.8K. In fact, it’s almost impossible to find sensible parameters to back up the “skating works 
by pressure melting” story. 

A salutary tale, all too often when a story is to good to check, even scientists don’t check. The 
concept is fine, but sometimes you just have to do the maths to show the application is wrong. 

And how does skating work? There’s probably friction and surface chemistry involved. The Fourth 
Law will explain it all... 
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10 Chemical Potential 


1. Chemical Potential: Nature’s boundary condition 

For a single component, the chemical potential is the specific Gibbs free energy, so start by con¬ 
sidering changes in that: 

dg = vdP — sdT 

Considering a constant temperature process, we have dT = 0 and 


dg 


dg m vdP = 



Doing the integral from reference state labelled “0”, we find: 


M - Mo = RT\n(P/P 0 ) 


Since the atmosphere is in contact with the ocean, the chemical potential of Carbon Dioxide is the 
same in each. Hence: 


rearranging 


,,atm ocean dt m f^atm /^ocean\ 

Mo - Mo = K1 in (Pi /Pi ) 

„ocean _ „atm // ,,atm ,,ocean\ /nm 

Pi - Pi expHMo - g 0 )/nl) 


with atmospheric pressure at 10 5 Pa, the partial pressure of CO 2 is 40Pa. Taking RT = 2500 Jmol 1 , 
and the given values for the chemical potential, we find that: 


p ocean = 40 exp(-8000/2500) = 1.4Pa 

Assuming that partial pressure is simply proportional to concentration, this implies that the con¬ 
centration of CO 2 in seawater is 0.0014% (14 parts per million). 

The difference in go is due mainly to potential energy: CO 2 in water disrupts the bonding network, 
increasing the potential energy of the solution by an amount proportional to the concentration of 
CO 2 . It is proportional because at 14ppm, the CO 2 molecules are far apart and so do not interact. 

2. Irn Bru 

Assuming that the chemical potentials come into equilibrium, both gas and drink have the same 
CO 2 chemical potential. We can take this as the ideal gas value: 

g = go + RT ln(P /Pq) = go + 5 RT 

when the cap is removed, the CO 2 partial pressure in the gas drops to its atmospheric value 
0.0004 x P 0 , so now in the gas 

g = g 0 + PTln(0.0004) = /i 0 — 7.824PT 

The release of pressure is done quickly, so the CO 2 in solution is the same as before. Gas and 
solution are out of equilibrium, so there is a significant difference in chemical potential between the 
two regions. CO 2 will fizz out of solution when shaken. A first approximation would be to treat 
the large volume of liquid as an infinite reservoir, so that the CO 2 would reach a partial pressure 
of 5atm (making 6 atm in all). This amount of CO 2 in gas would be 

_ PV _ 5 x 10 5 x 0.05 x 10 ~ 3 
n ~ H/T ~ 8.314 x 300 


which is 0.01 moles 

We were told that originally there were 0.155 moles of CO 2 , so the approximation that most of 
the CO 2 remained in solution is a good one: the chemical potential in the drink will, in fact, have 
dropped by about 6 %. 
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Repeated iterations will extract progressively less CO2, the concentration decreasing by a geometric 
series. We could calculate this, but since there is no numerical definition for “flat” it is sufficient 
to use the 6% value to estimate that after 10-20 iterations, most of the CO2 will be lost. 

The -7.824RT was never used. This is equivalent to neglecting the additional atmospheric CO2 
added to the bottle on each iteration. Notice also that if there were no CO2 in the atmosphere, 
the chemical potential would be minus infinity : the gas is infinitely attractive to the first CO2 
molecules to arrive, or as Aristotle put it “Nature abhors a vacuum”. This infinity does not mean 
that the Gibbs free energy is infinite, because the free energy of any system is the chemical potential 
times the number of particles in the system. This product goes as nlnn, which goes to zero as 
n —> 0. 

3. Chemical Potential Change in Mixing 

Once the containers are connected, we are out of equilibrium, and expect the gases to mix. Assume 
that both argon and krypton can be treated as ideal gases, which do not interact. 

From the First Law, no work is done and no heat enters the system, so the internal energy cannot 
change. Since the internal energy of the ideal gas depends only on temperature, it is initially the 
same for each component. At equilibrium, the system has a single temperature, so we can deduce 
that the temperature cannot change. 

The volume doubles for each gas, at constant temperature, so since we have an ideal gas, the partial 
pressure halves. However, since the two gases intermingle, the total pressure is unaffected. 

The entropy change is given by 


ds = °^-dT + —dv 


where v is the specific volume, 
argon: 


This is a quantity which changes by a factor of 2, hence for the 


dS = R 


ln2V 

V 


= R In 2 


and the entropy of the krypton changes similarly. So the total entropy change of the universe is 
2Rln2. The same value can be calculated from microscopic considerations, that each atom could 
be one the two types, so for 2Na molecules: 


S = k B \nW = k B \n2 2NA = 2Rln2 


where Avogadro’s number is Na = R/k B . Note that the thermodynamic derivation uses only 
macroscopic quantities, so is independent of the existence of atoms. The Gibbs free energy (chemical 
potential) for each species changes by: 


A gi = A Ui — A (Tsi) + A (piVi) = 0 — RT In 2 + 0 

Notice that we use the partial pressure, which is halved as the volume doubles. No work is done 
in the mixing, so we should expect the p^i term to make zero contribution to the change in free 
energy. 

If the initial volumes were different, then the number of moles of the two gases would be: 2Nk /(W4 + 
Nk) and 2 Na/{Na + ); the entropy changes of the two gases would be different, but still given 
by the same equation: 


2 N k R (V a + Vk) ZNaR V a + V K 

{N a + Nk ) n Vk ( N A + N K ) n Va 

In terms of mole fractions, this would be 

dS = —2 R{xk In xk + %a 1 n^) 

Notice that if there is only one gas {xk = 1; xa = 0) there is no entropy of mixing, dS = 0. 
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4. Regular solution and solubility limits 

This expression for the potential energy is proportional to the densities of each gas, so represents 
the average number of interactions between different types of atoms. This is similar to the a/v 2 
term in the van der Waals equation, but of opposite sign because they are repulsive. 

In the isothermal process kinetic energy (temperature) is conserved. The potential energy rose by 
z = ZxaXk , so the internal energy increases by the same amount, 


Ag = Au + TAs + APv 

The change in entropy and Pv is the same as in the previous example (except the total is now 1 
mole, so xa+ %b = 1): so we can write: 


A g = ZxaXb — RT(xb In xb + xa In xa) 


The most stable structure is always the one which minimises A g, i.e. 


dAg 

dxA 


= 0 



For high temperatures, relative to the repulsion Z, the lowest Gibbs free energy can be read directly 
from the graph. For lower temperatures, there are two minima, and the lowest Gibbs free energy 
comes from a phase separation into A-rich and B-rich regions, the amount of each region being 
determined by the overall concentration via the lever rule. 

Differentiating gives: 


= Z( 1 - 2x a ) + RT ln(au) - RT ln(l - x A ) = 0 
ax a 

which (via symbolic algebra package for Z=3RT) has solutions at x=0.07, 0.5 and 0.93. The first 
and third are minima. 

An approximate solution to the equation can be obtained by assuming that for small xa' In (xa) — 
ln(l — xa) ~ In (xa)\ (1 — 2 xa) ~ 1 • Now 

Z( 1 — 2 xa) + RT\h{xa) = 0 
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or 


xa = exp(— Z/RT) 


This is always positive, so no matter how strongly the fluids repel, some solubility is always ther¬ 
modynamically favoured. Notice also that solubility always increases with increasing temperature. 
There are limiting cases at xa = 0,1. These are not minima of the function but correspond to 
physical limits on concentration xa • For an unphysical concentration, the equation for g gives a 
complex value. 

Aside: xa = exp {—Z/RT) looks the same as the Boltzmann probability in statistical physics for 
finding a microstate with excess energy Z. But notice that the derivation from thermodynamics 
considered only macroscopic quantities. 

The chemical potential in a mixture is given by the value of g. Substituting in the values at the 
turning points gives: 

gA = -0.05 8RT, +0.057 RT, -0.05 8RT 

The stable state is the one with minimum chemical potential. 

By symmetry, the chemical potential of B will be the same as for A. Up to x=0.07 this is a mixture 
of the two substances. For concentrations between 0.07 and 0.93, the lowest Gibbs free energy is 
obtained with a two-phase mixture of an A-rich and a B-rich fluid. Above x=0.93 a single phase 
is obtained. x=0.07 is the solubility limit. 
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5. Supercool 

Firstly, we calculate the • Assume that we have an ideal solution at atmospheric pressure, 

then the partial pressure of water changes from P 0 to Po(l — X), and using the expression for an 
ideal solution, the chemical potential is reduced by 

/i(X, T) - 0, T) = RT ln(l - X) 


So that 


dfiL \ _ —RT 

~dX ) t ~ l-X 


The general expression for a change in chemical potential at constant pressure is: 



dT + 
p,x 



dX 

P,T 


Since dg = dg 


—sdT + udP, we can write 



p,x 


—s 


and therefore 

d(i = —sdT — RTdX 

On adding the salt, we must have the same change for liquid and ice dfiL = d/i/. Since dX = 0 for 
the ice, equating the chemical potentials gives: 

dfij = —RTdX — sidT = d/ij = —sjdT 


The entropy difference is simply related to the latent heat by (si — s/) = Z/T, so we now have 

dT _ RT 2 _ 8.3 x 273 2 _ 

dX ~ T ~ 18000 

so with a change AX = 0.1 we have a 3K drop in temperature. 

In fact the temperature drop is bigger than this - the chemical bonding between water and salt 
lowers its chemical potential by more than the entropy of mixing. Notice that this drop in temper¬ 
ature is not due to a loss of energy (e.g. in an endothermic reaction). Some of the ice has melted 
and the latent heat absorbed in this process was extracted from the liquid. 


6. Simplified Osmosis 

Consider the two systems, inside and outside the cell. Since water molecules can pass through the 
cell wall, the chemical potential of water must be the same on both sides. Protein molecules cannot 
pass through the semipermeable membrane, so their chemical potential may be different. 

Treating water and protein as ideal gases, the partial pressure of the water inside the cell is the 
same as the water pressure outside, latm. So the difference in pressure is just the partial pressure 
of the protein, which is 2% that of the water. Since the creature lives at the surface, the external 
pressure is latm, so the excess pressure is 0.02atm. 

If we consider a hemisphere, then the net force due to excess pressure is balanced by the stress in 
the membrane. 


A P x 7rr 2 = a x 2tt rt 

Rearranging to get stress, with diameter 2r = 10 -5 m and thickness 10 -8 m 

cr = 0.02 x 10 5 x 10“ 5 / 1( r 8 = 2 x 10 9 Pa 
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The Essential Mathematics: 

1. Partial differentials applied to an ideal gas 

The volume in an ideal gas is a function of two variables, e.g. V(T,P). To fully evaluate the 
derivative, we need to know what quantity is being help constant. An isobaric expansion is different 
from an adiabatic one. 

The ideal gas equation is PV = nRT , so the thermal expansion coefficient at constant P is 

1 f dV \ _ 1 nR _ 1 
V [dfjp ~ Vi 5 " ~ T 

For the adiabatic expansion, we want to replace the explicit P by a constant, say K = PT 7 , so 
the ideal gas law is 

KV 1-7 = nRT 

and the adiabatic thermal expansion coefficient is: 

1 fdV\ _ 1 /niA (1/(1-7)) d(r 1 /( 1 -T')) _ 1 

V\df) K ~V\YJ dT _ (I -7 )T 

Which is actually negative! The two “thermal expansions” are different because they refer to 
different processes: in the first, familiar, case we describing an expansion due to heating without 
changing the pressure. The adiabatic case is a peculiar situation where we increase the temperature 
without supplying any heat. For an ideal gas, that means increasing the pressure and thus reducing 
to volume. Since 7 > 1 we see that (f^) P is positive, but K is negative for an ideal gas. 

2. Exact differentials 

Given X = PV 3 + aT, the constant a must have units of Jm 6 /K. For an ideal gas PV = nRT so 
X = T(nRV 2 + a) 


and 

(w)r = 2nRVT 

Wv(wl = 2nEV 

The second derivative is the same, regardless of order of differentiation, hence X may be a state 
variable, although not a very useful one! 

For Work, we can consider 

fdW\ 

(ap)v = ° 

since no work is done at constant volume. However 

_d_ / dW\ _dP 
dPv\~dV)p ~ dPv ~ 

Thus the second differential d pJ^ Vp does depend on order of integrations, and so work is not a 
state function. 
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